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EXPERIMENTAL  OBSERVATION  OF  SMALL-SCALE 
THERMAL-BLOOMING  INSTABILITY 

Zhang  Kai ,  Shi  Jiang jun,  Chen  Junyan,  Qiu  Fumin r 
Yue  Zhengpu,  and  Yang  Chenglong 


Southwest  Institute  of  Fluid  Physics 
P.  O.  Box  523,  Chengdu  610003 


ABSTRACT  The  experimental  simulation  results  of  instability  of  Stimulated  Thermal  Rayleigh 
Scattering  { STRS  )are  reported  in  the  paper,  which  include  the  studies  on  the  process  of  instabili¬ 
ty  of  the  sinusoidal  perturbation  source  propagating  in  the  thermal  -  blooming  medium  and  the 
experimental  studies  on  the  increase  of  the  small-scale  thermal  -  blooming  instability  under 
different  perturbation  frequencies  and  different  light  power,  and  thus  proves  that  the 
small-scale  instability  caused  by  thermal  -  blooming  increases  exponentially  as  laser  propagation 
in  the  absorbing  medium. 

KLY  WORDS  laser  beam  propagation  in  the  atmosphere,  thermal  -  blooming,  small-scale 
thermal  -  blooming  instability. 


I.  Introduction 

In  1986,  Jan  Hermann  was  the  first  to  discover  rapid  growth 
in  small-scale  thermal -blooming  instability  during  laser 
transmission  in  computer  simulation.  Thereafter,  theoretical  and 
experimental  studies  on  small-scale  instability  became  the  focus 
of  thermal-blooming  research  [1-7].  Small-scale  instability 
includes  stimulated  Rayleigh— scattering  (STRS)  instability, 
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stimulated  thermal— Brillouin  scattering  (STBS)  instability 
phase -complement  ins t a^"^  i  u  +-T7-  /dp t ^  Tp  f tn^-v-t/-^  r\ -f  iqqq  +- 

1990 ,  Lincoln  Laboratories,  the  Northeast  Research  Conference, 
and  the  Livermore  National  Laboraf nv*T7  vonnrfori  r3n-i  hi 

growing  phenomena  of  PCI  and  STR.S  instabili tie°  pao^aw-h-rc 

generally  consider  that  the  probability  of  small-scale 
instability  may  be  the  final  limiting  condition  of  threshold 
value  of  transmission  intensity  of  high-powered  lasers  in  the 
atmosphere  with  respect  to  the  reflective  surface  of  high-powered 
lasers  with  large  aperture  and  short  wavelength. 

This  paper  reports  that  the  authors  utilized  an  argon  ion 
laser  (single  wavelength  with  maximum  output  power  1.5W)  to 
simulate  thermal-blooming  cells  by  a  carbon  tetrachloride  and 
iodine  solution.  The  experimental  results  of  rapidly  growing 
perturbation  bands  of  intermediate  and  small  spacing  scales  in 
t  h.  e  t  h  e  r  m  a  1  ~~  jb  1  o  o  mi  1  n.  oj  f  o  r  m  a  1 1  o  m  ^  v*  <■— \  n  o  c  c  3  v*  q  ^  o  c?  o  t-s  ^  ^ 
Experimentally,  the  existence  of  small-scale  instability  was 
proven.  Moreover,  some  experimental  studies  were  con^ 1 1  n ^ Qd  An 
the  growing  role  of  small-scale  instability  at  different  light 
power  and  different  perturbation-spacing  frequencies. 


II.  Principle  of,  and  Optical  Path,  in  the  Experimental  System 
Fig.  1  shows  the  principle  of  the  experimental  system.  In 
the  experiment,  the  wavelength  of  the  argon  ion  laser  is 
0 . 488mum,  operating  at  a  single  wavelength.  The  output  power  is 
continuously  adjustable.  When  the  output  power  is  relatively 
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Fig.  I  Layout  of  the  STRS  experiment 

argon -ion  laser  ;  2.  shutter  ;  3.  4.  17.  18.  beam  expension  mirrors  ;  5.  8.  reflectors 
6.  flitter  ;  7.  screen  1  ;  9.  pertubation  source  :  10.  thermal  -  blooming  cell 
11.16.  beam  spliter  ;  13.  screen  2  :  14.  CCD  camera  ;  19.  He-  Ne  laser 


small,  the  output  mode  is  single  fundamental  mode  (TEMsubOO) 
The  diameter  of  the  thermal  blooming  cell  is  9cm;  its  length  is 


cn  i  on. 


The  thermal  blooming  medium  is  carbon  tetrachloride.  The 
absorption  coefficient  of  a  light  beam  passing  through  the 
thermal  blooming  cell  can  be  controlled  by  adding  small  amounts 
of  iodine.  The  result  of  growing  small-scale  instability  of 
thermal  blooming  was  recorded  with  a  high-speed  camera  and  a  CCD 
solid-state  camera.  After  being  gathered  by  the  image  collecting 
system,  the  experimental  data  were  processed  on  a  computer.  The 
fastest  sampling  frequency  of  the  high-speed  camera  is  144  frames 
per  second.  The  CCD  real-time  image  collecting  system  collects 
32  frames  of  128x128  images  (or  8  frames  of  256x256  images)  at  a 
speed  of  16.6  frames  per  second.  The  sampling  time  of  the  CCD 
camera  is  controlled  with  an  electronic  shutter;  the  shutter 
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light  integration  time  is  1/10 ,000-th  of  a  second.  The  beam 
expansion  lens  complex  adopts  the  composition  of  positive  an 
n  e  g  a  t  a  v  e  1  e  n  s  e  s  n  o  (j  0  t'*  "f-  a  ^  ri  i  n  c  ^  ^  ^  ^ 

the  straightening  properties  c 
mall-scale  perturbation  source 
ass  plates  with  a  certain  ang! 

The  perturbation  source  is  si 
o w  of  the  thermal -bl  ooming  cs] 
ough  the  perturbation  source  e 
f  t  h  e  g  lass  plates,  the  V)  0  03  i 
laser  beam,  generating  perturbation  bands  with  a  distribution  of 
a  certain  sinusoidal  intensity.  The  spacing  period  of  ^he 
can  be  adjusted  by  the  ang 1 e  me luded  between  two  surfaces  of  th 
1 1  glass  plates ,  and  the  placement  of  the  band  in  the 
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optical  path. 

Experimental  light  path:  after  the  output  of  the  argon  ion 
laser,  the  motorized  shutter  controls  the  light  exposure 
synchronization  and  the  exposure  time,  for  synchronizing  with  the 
CCD  camera  and  the  image  collection  system.  After  beam  expansion 
of  the  laser  beam  through  the  lens  complex,  it  becomes  a 
straightened  beam.  After  passing  through  the  perturbation 
source,  the  beam  passes  through  a  thermal— blooming  simulation 
cell.  In  the  exit  optical  path  an  imaging  lens  ic  ^  1  a ^ -1- ^ 
converge  the  light  intensity  images  at  the  exit  of  the  thermal 
blooming  cell,  to  form  an  image  on  screen  1,  then  the  image  is 
recorded  with  the  high-speed  camera  or  CCD  camera  system.  In  the 
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experiment,  another  beam  of  the  He-Ne  laser  x s  ejcjpandsd  in  tine 
reverse  direction  and  passes  through  the  thermal  blooming  cell. 

A  filter  plate  isolates  the  He-Ne  laser  in  the  main  optical  path 
to  be  projected  onto  screen  2.  The  beam  interferes  with  another 
beam  of  reference  light  for  interference  measurements  for 
observation  of  the  stability  of  the  carbon  tetrachloride  solution 
in  the  thermal-blooming  cell.  The  observation  coefficient  of  the 
thermal  blooming  cell  with  respect  to  the  argon  ion  laser  beam  is 
obtained  by  measuring  the  power  at  the  entrance  and  exit  of  the 
thermal  blooming  cell  with  a  laser  wattmeter. 


III.  Experimental  Parameters  and  Experimental  Results 
1.  Typical  experimental  conditions: 
laser  wavelength  A=0 . 4 8 Smicrometer ; 
absorption  coefficient  a=0.05cm'*; 
laser  output  power  P=366mW? 

perturbation  period  (of  small-scale  perturbation  source) 
d=Q . 0864cm ; 

after  laser  diffusion,  light-spot  diameter  (l/e“), 

D=0 .68cm; 

thermal-blooming  cell  length  L=60cm; 
initial  perturbation  intensity  A///%±0.0 8 


2.  Calculation  of  parameters 
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Fennell  number  (of  the  entire  beam)  A /  —  ~  *  —  83.8  ,  ^  1.48 
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thermal-blooming  growth  rate  — - —  =  12.76s'1 

for  small-scale  perturbation  the  characteristic  expansion 

time  Tc=  (DK2)~'=* 0.234s  ;  D=  =8.1  x  10'4cnr/s 
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3  #  Experimental  r ssul t e 
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From 

Fig .  2  f  we  c< 

Fig.  2  Spatial  distribution  of  argon  -  ion  laser  spot 
intensity  at  the  thermal- blooming  cell  exit 


see  that  the  light  intensity  distribution  is  not  the  ideal 
Gaussian  distribution  due  to  the  effect  of  fluctuations  in  lase 
scattered  spots.  The  contrast  of  the  prearranged  perturbation 
bands  is  smaller  than  the  estimate  based  on  theoretical 
calculations.  Fig.  3a  shows  the  light  spot  intensity  lines 
passing  through  the  center  line  and  perpendicular  to  the 
perturbation  bands  when  t=0 .  Fig.  3b  through  3f  show, 
respectively,  the  intensity  lines  beginning  at  thermal  blooming 
then  at  time  intervals  60,  120,  180,  240,  and  300ms.  At  the 
beginning  of  thermal  blooming,  due  to  the  flaw  points  of  the 
laser  scattered  spots  and  the  optical  system,  the  prearranged 
perturbation  bands  are  not  distinct.  With  prolonging  of  therma 
blooming  time,  we  can  see  that  the  contrast  of  the  perturbation 
bands  clearly  increases.  This  proves  that  during  transmission 
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Fig.  3  Intensity  distribution  at  center  line  at  time  0.  60.  120.  ISO.  240,  500ms 


3  lcissir  beam  in  3  medium,  the  small-scale  =-i- 

certain  spacing  frequency  rapidly  grow  with  the  generation  of 
thermal  blooming. 

Fig.  4a  through  4f  indicate,  respectively,  the  Fourier 
transform  frequency  spectra  of  light  intensity  distributions  at 
time  0,  and  after  the  beginning  of  thermal  blooming,  at  time 


intervals  of  60,  120,  180,  240,  and  300ms.  The  corresponding 
coordinate  point  of  the  spacing  frequency  spectrum  of  the 
prearranged  small-scale  perturbation  is  1.3mm"'.  From  the 
corresponding  frequency  spectral  components,  we  can  see  the 
obvious  growth  in  small-scale  perturbation  bands. 

Fig.  5  shows  the  growth  process  with  time  in  the  spacing 
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frequency  spectral  components  of  small-scale  perturbation. 
Figs.  6a  through  6f  show  a  series  of  pictures  of  small-scale 


+-  H 


hermal  blooming  instability  growth  at  the 


corresponding 


intervals . 

Moreover,  the  authors  conducted  increasing  experiments  on 
small-scale  perturbation  of  different  laser  power  and  the  width 
of  different  prearranged  perturbation  bands.  Table  1  shows  the 
experimental  conditions  and  the  results.  By  comparing  the 
experimental  results  of  several  experimental  conditions,  we 


obtained  the  following: 

1.  For  the  same  laser  power,  for  different  small-scale 
perturbation  periods,  the  characteristic  time  of  diffusion  is 
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Table  I  Comparison  of  the  small-scale  instability  increasing 
under  differeot  experimental  condition 


E 

T,-(D  ■  K1)-'/  s 

'■-JJh 

small  ~  scale 

instability 

No.  0 

366 

0.0864 

0.234 

7.97  x  I0'> 

obvious 

No.  1 

184 

0.154 

0.742 

1.12  x  10'1 

obvious 

No.  2 

594 

0.069 

0.149 

6.25  x  JO*’ 

not  obvious 

No,  3 

73 

0.154 

0.742 

1.78  x  |0'J 

not  obvious 

No. 4 

594 

0.108 

0.365 

6.25  x  i<r> 

obvious 

Fig. 5  The  small-scale  perturbation  spatial  spectrum  amplitude  vs  lime 

(a  )  r  =  Oms  ;  (b),  =  60ms;  (c)r-  120ms  ;  (d)r-  ISOmr;  <e)/-240ms  ;  (f),  =  300ms 


i  n 


Hg.6  A  series  of  spot  pictures  taken  at  different  time  (a)  —  (f)  f~0,  60,  120,  180,  240,  300ms 


short  for  small  perturbation  periods  (No.  2) .  The  medium 
diffusion  effect  plays  a  major  role,  exhibiting  the  absence  of 
gam  for  the  small-scale.  For  large  perturbation  periods 


(No.  4) ,  the  diffusion  characteristic  time  is  long.  Now  th 


e  STRS 


instability  grows ,  playing  a  major  role.  This  exhibit 


c  p  r  p.pv 


growth  in  small-scale  prearranged  perturbation. 

2.  For  the  same  small-scale  perturbation  period,  with 
different  laser  power ,  in  experiments  (No.  1)  at  large  laser 
power ,  the  small-scale  prearranged  perturbation  bands  clearly 
grow.  Conversely  (No.  3)  for  low  light  power,  in  this  small- 
scale  perturbation  spacing  period,  the  diffusion  effect  of 
refractivity  fluctuations  of  the  medium  plays  a  major  role.  This 
shows  that  small-scale  perturbation  does  not  have  the  phenomenon 
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3.  The  experiment  proves  that  the  small— seal 

perturbation  period  has  a  very  great  effect  on  the  rapid  growth 
in  small-scale  thermal-blooming  instability.  This  is  consistent 
with  the  indicator  terms  in  F^^  f q  ^  m  o  nmnArf-innai  +- r\ 

“  *  *  4  '■  -L^  Vi.  U4.  WAAU  J.  ^ 

the  square  of  K. 

4.  In  the  thermal  blooming  process ,  whether  or  not  the 
contrast  of  small  —  scale  perturbation  band  intsnsi^”"^  00  i o  f v» q 
competitive  result  of  the  diffusion  effect  that  small-scale 
perturbation  rapid  growth  due  to  light  transmission  thermal 
blooming,  and  the  diffusion  effect  of  the  refractivity 
fluctuations  in  the  medium. 


IV.  Conclusions 

In  the  small-scale  instability  experimental  studies  on  open- 
ring  thermal  blooming/  with  an  argon  ion  laser  passing  through  a 
simulated  thermal-blooming  cell  with  carbon  tetrachloride  and 
iodine  solution,  the  authors  observed  the  growth  of  sinusoidal 
intensity  perturbation  bands  in  different  experimental 
conditions.  This  illustrates  that  the  rapid  growth  of  small- 
scale  perturbations  is  mainly  determined  by  the  growing  rate  of 
thermal  blooming,  and  the  magnitude  of  the  spacing  period  of  the 
small  spacing  seal© 

The  first  draft  was  received  on  February  25,  1993;  the  final 

_ v e d  for  publication  on  June  1  1993 
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ABSTRACT  Through  actual  atmosphere,  we  have  measured  atmospheric  transmittance  for 
1.315/im  Iodine  laser.  With  medium  visibility,  relative  humidity  64%  and  2lC  .  average  value  of 
atmospheric  transmittance  is  0.78.  This  result  is  fairly  identical  with  computation  value 
obtained  by  existing  data  and  computer  program. 

KEY  WORDS  iodir.e- laser,  atmospheric  propagation,  atmospheric  tranmittance. 


I .  Introduction 

The  1 . 3 1 5— micrometer  iodine  laser  located  at  the  window 
region  of  the  atmosphere  evoked  great  interest,  since  it  has 
)us  potential  as  a  high-energy  laser.  Many  specialists 
search  on  the  laser  adsorption  in  1  sbor a ^ ^ o  -f 
water  vapor,  carbon  dioxide,  methane,  and  other  gaseous  molecules 


AV'rrtm 


cl  n  ri  f  csrl 


1  A 


Fig.  I  The  schemalic  diagram  of  experiment  system 


II.  General  Experimental  State 

Fig.  1  is  the  schematic  diagram  of  the  entire  experimental 
system.  The  photclytic  iodine  laser  device  v?as  developed  and 


provided  by  the  Dalian  Institute  of  Chemistry  and  Physics  of  the 


Chinese  Academy  of  Sciences.  The  pulse  width  of  the  iaser  wa0 
lOOmicroseconeds ,  since  the  outputted  laser  beam  was  of  ring- 

shaped  low — order  mode.  The  exterior  diameter  ^ e  -h  v*  ^  _ ch apori 

light  beam  was  18mm,  and  the  outputted  laser  energy  was  8 Om-J .  In 
addition  to  its  use  as  light  guidance,  the  He-Ne  laser  was  also 
used  as  a  light  source  to  measure  the  atmospheric  transmittance 
of  the  0 . 6328~  micro  meter  laser.  The  atmospheric  transmitt  an^0  ^  -f 
t  h  e  w  a  v  e  1  e  n  g  t  h  w  a  s  u  sed  to  esti  m  ate  the  ^  "i-  t*>  /—*  n  onf  o*  ■f 

the  atmospheric  aerosol  particles  in  the  total  experimental  liaht 
path.  Ho w ever,  lasers  at  both  wavelenoths  were  transmitted 
through  a  Cassegrain  telescope  with  an  aperture  3mm  OD .  At  the 
receiving  end,  there  was  an  aluminum  coated  concave  mirror  ( m ~ 1 m ) 
with  OD  290mm.  After  convergence  of  the  received  laser  beam,  a 
laser  wattmeter  was  used  to  conduct  detection  of  1 . IS  — micrometer 
laser  energy.  The  0.S328  —  micrometer  laser  signal  t.tsc? 
with  a  photomultiplier  tube.  To  obtain  the  real-time  data  of 
atmospheric  transmittance,  synchronous  monitoring  was  conducted 
at  the  transmission  terminal  on  the  laser  transmission  energy  or 
power  of  both  wavelengths  by  using  the  laser  wattmeter  and  a 

after  measurement  experiments  each  day,  comparative  calibration 
of  the  system  was  made  on  the  instruments  at  both  terminals.  By 
definition,  the  calibrated  coefficient  of  the  system  multiplied 
by  the  ratio  of  the  measured  values  between  the  receiving  end  and 
the  transmission  end  is  the  actual  atmospheric  transmittance  of 

the  franem-i  con  nn  1  -i  +-  r-\  rs  +-  Vi 
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The  transmission  distance  bet 
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the  transmission  end  and 
ie  Q  .  46km .  Although  it  v/as 


mid — autumn  ygt  mostly  cloudy  during  the  exp)erimental  period, 

t*t  "i  -H  Vi  f  rammn  f-  v  *i  nr  rr  1  tt  r»r  €  r^rrrr'tr  ^a^re  rP'Hcj  p  *1  a  f  i  t;  o  V-itim-i  r?  -i  f  rjae 

Vi  X  VU  -X.  v_^  V—F  X  X  U  VA  X-  X  Lj  ZJ  -L  J  i.  ^  ^  ^  y  J  ^  V*  Jf  •  J.  X  A  X  A.  SAiU  -i-  VA  x  W  j  VI  V-A. 

between  SO  and  80%.  There  were  rarely  sunny  days.  To  ensure  the 
correct  processing  of  the  experimental  data,  the  authors  made 
observations  in  the  entire  experimental  period  as  well  as  to  the 
conventional  meteorological  parameters  (such  as  temperature, 
relative  humidity,  atmospheric  pressure,  and  visibility  in  the 
atmosphere) . 


III.  Data  Processing 

During  laser  transmission  in  the  atmosphere,  before  the 
emergence  of  the  nonlinear  effect,  the  transmittance  can  be 
described  with 


T  =  exp[  - 


/i  ( r)dr ] 


(1  ) 


which  is  Lambert's  law.  In  the  optical  path  of  horizontal 
transmission,  the  atmospheric  medium  in  the  optical  path  can  be 
considered  as  a  homogeneous  medium.  Therefore, 


7=exp[  -fiL] 


(2  ) 


In  Eq.  (2),  L  is  the  transmission  distance;  and  p  is  the 
attenuation  coefficient  of  the  atmospheric  medium.  The 
coefficient  includes  the  absorption  coefficient  cr  and  the 
scattering  coefficient  be tam  induced  by  absorption  and  scattering 
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^jr  the  gas  molecules  in  the  s tirosphsre  r  es  well  s. s  the 
coefficient  due  to  absorption  and  s cat ter i no  of  the  aerosol 
particles  in  the  atmosphere.  That  is,  the  attenuation 


nriof  f  ■ 


c  i  e  n  t 


(3) 


In  actual  atmospheric  transmission  of  the  1 . 315-micrometer 
iodine  laser,  besides  the  aerosol  particle  scattering,  absorption 
by  water  vapor  molecules  in  the  atmosphere  is  the  main  factor 
sf  f  6ctiriy  transmittance  •  There  are  other  atmospheric 
ps^ticles  (such  as  carbon  dioxide,  methane,  and  ammonia)  also 
making  a  certain  contribution  to  the  atmospheric  transmission 
attenuation  of  the  iodine  laser  wavelength.  However,  this  is  a 
small  quantity  compared  to  the  absorption  of  water  vapc*- 
molecules.  Therefore,  in  this  article,  we  do  not  consider,  for 
the  time  being,  the  attenuation  function  for  a  1.315-micrometer 
laser  absorbed  and  generated  by  these  gas  molecules.  During  data 
processing,  the  functions  due  to  these  gas  molecules  are  also 
considered  as  the  absorption  coefficient  of  atmospheric 
molecules.  At  the  1.315-micrometer  wavelength,  the  molecule 
scattering  coefficient  is  relatively  small  and  can  be  neglected 
In  actual  atmospheric  transmission  of  an  iodine  laser,  the 
mam  function  affecting  the  transmittance  is  scattering  of  the 
atmospheric  aerosol  particles,  and  absorption  by  the  atmospheric 


1 0  1  Pm!  Q  C 


In  the  authors '  experiments ,  the  atmospheric 


transmittance  of  the  iodine  laser  obtained  can  be  wr i t t* 


;n  as 


1  Q 


r,=exp[-(o.+  ^)L] 


(4) 


Here  we  neglect  the  molecule  scattering  coefficient  {3*  at  the 
1 • 315 — micro meter  wavelength 

With  respect  to  the  radiation  transmission,  by  using  the 
scattering  of  atmospheric  aerosol  particles  as  the  main  light 
extinction  factor,  generally  the  following  empirical  formula  i i 
used  to  estimate  the  extinction  coefficient. 


H  =  0.912/ VU)(0.S5/X)' 


(5) 


In  the  equation,  V»  is  atmospheric  visibility  (km)  ;  X  is  the 
radiation  wavelength  (micrometers) ;  and  q  is  the  visibility- 
wavelength  revision  factor. 


0.585  V"' 


6km  <  Ky<20km 
!■  KJ/>20km 


(6) 


The  effectiveness  of  Eq.  (5)  has  been  proven  by  a  series  of 
experiments  by  the  authors  [4] .  The  attenuation  of  an  He-Ne 
laser  for  A=0 . 6328micrometers  is  the  result  of  the  scattering 
function  of  the  atmospheric  aerosol  particles.  However,  the 
absorption  and  other  functions  of  the  atmospheric  gaseous 
molecules  can  be  neglected.  Therefore,  we  can  use  the 
atmospheric  transmittance  of  the  0.6328-micrometer  laser,  and 
estimate  the  extinction  coefficient  for  scattering  by  the 
atmospheric  aerosol  particles  with  Eq.  (5).  If  in  the 
experiments  the  atmospheric  transmittance  is  Ty  for  the  0.6328- 
micrometer  laser,  and  the  scattering  coefficient  for  the  aerosol 
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particles  is  {3n-n ,  then  we 


:an  obtain  the  following"  from  Egs 


/  o  \ 

\  ^  i 


(  5  ^  an^ 


( a ) 

\  w  / 


P,n  =  ~  (In Tw)/L  =  (3.912/  )( 0.55/ 0.6328  V 


(7) 


Py  ciiK 


c?  /dm  o 


s  t  i  tut  mg  A,— 1 . 31 Bmicrometers  into  Eg.  (5)  r  with 
manipulation  we  can  obtain  the  scattering  coefficient  of  the 
iodine  laser  for  the  aerosol  particl 


O  O  O  c? 


Pip~  ft ,11  (0.6328/  1.315  )*«=  ”[ (In Tw )/£] (0.4812  )f 


(8 


Substitute  Eg.  (S)  in  Eg.  (4)  and  we  he 


T,  =  exp[  -  ( a„+  p,F )  L]  =  exp{  -  ( a.-  0.48 1 V  ( In T  „ ) /  L)  L  }  ( 9  ) 


T]  is  the  atmospheric  transmittance  of  an  iodine  laser  in  actual 
measurement.  Thus  w e  can  obtain  the  absorption  coefficient  of 
the  iodine  laser  due  to  the  atmospheric  molecules  as 


=  [  (0.48 12  )*  InT^-lnT,]/  L 


(10) 


During  the  horizontal  atmospheric  transmission  ( L~0 . 46km )  of  the 
iodine  laser,  the  atmospheric  transmittance  T  (km-*)  can  be 
written  as 


r=exp{-[(-lnr„)/0.46]}  (11) 

Wh en  we  ^ 
absorption 
be  used  to 

T  '=  exp{  [0.481 2*  In  T,  —  In  Tlt]  /  0.46 }  (12) 


nsider  the  transmittance  T 1  (km  *)  durin^  ths 


of  atmospheric  molecule0  +- q  f  i  ami  nr 
describe  the  transmittance  * 


eguation  can 


on 


In  the  equation,  the  value  of  q  is  determined  by  (6) 


anno 


IV.  Results  and  Discussion 

Table  1  lists  some  of  the  typical  data,  including  the 
atmospheric  transmittances  (Tp  and  Ti )  measured  at  the  dist; 
of  460m  for  an  He  - IMe  laser  and  an  iodine  laser.  Derived  from 
Eqs .  (11)  and  (12),  the  atmospheric  transmittance  of  the  iodine 

laser  per  km  (T)  and  the  per-km  transmittance  (T')  in  the 
situation  of  absorption  by  molecules,  with  deduction  of  the 
aerosol  extinction  coefficient,  we  obtain  the  molecule  absorption 
coefficient  a„.  The  first  three  columns  are  the  visibility  V« , 
relative  humidity  (H) ,  and  temperature  (t) . 

From  the  data  in  Table  1,  we  can  see  the  following:  during 


Table  1  The  typical  measurement  results  of  atmospheric  transmittance 
for  iodine  laser  (1991.10) 


Date 

91  -  10 

Time 

VJ  km 

HI  % 

t!  t 

■ 

■ 

*m!  km"' 

6 

09  : 30 

7.0 

0.800 

0.878 

0.754 

0.909 

0.095 

09  : 30 

7.0 

0.800 

0.863 

0.727 

0.876 

0.132 

09  : 30 

7.0 

0.800 

0.869 

0.736 

0.888 

0.1  19 

z 

15  :  10 

7.0 

0.800 

0.850 

0.701 

0.846 

0.167 

15  :  15 

7.0 

'  0.800 

0.905 

0.806 

0.972 

0.029 

1 5  :  30 

7.0 

0.800 

0.888 

0.772 

0.932 

0.071 

15  : 40 

7.0 

0.800 

0.864 

0.729 

0.879 

0.129 

15  :  53 

7.0 

0.800 

0.840 

0.685 

0.827 

0.190 
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VJkm  H]  % 


1870 

0.903 

0.801 

0.93! 

0.857 

0.921 

0.836 

0.933 

0.860 

).870 

0.928 

0.850 

).870 

0.909 

0.813 

)  .8  70 

0.934 

0.862 

1.870 

0.894 

0.783 

1.870 

0.920 

0.834 

0.880  0.881 
0.880  0.885 

0.880  0.877 


w j_c  transmission  of  a  1 .  BIS-inicroiseter  iod^r^  i  ^ ^ ^ 

transmittance  varies  mainly  with  the  variaf  ■*  °  ^  r  af^cphc‘v'-!  ^ 

.^i^ili  ty  and  water  vapor  content.  The  average  meteorological 

conditions  are  VH=10km,  H=64%,  and  t=21cC.  For  the  iodine  lasei 

the  averaae  transmi  n  7  q  u  ™  ~  ~  n  *  u  _ 

^  wiiuiux  ^  vdii^v,  j.  yj .  /  o Km  ,  and  the  corresponding 

attenuation  coefficient  p=Q .  246km"' .  Among  the  latter  *-ha 

molecular  absorption  coefficient  o„=0 . 107km'1 ,  and  the  aerosol 

extinction  coefficient  j3n=0 . 139km-l 

To  verify  the  trustworthiness  of  the  above-mentioned 

results,  the  authors  made  computations  by  using  !*<=>  =  k-i^t 

resolving-power  spectral  lines  of  AFGL  in  the  United  States,  an 

the  LASER  program.  Table  2  shows  the  result*?  tv-,,,* 


thst  if  t h 6  two  types  (urban  type  and  rural  type) 


=)  £=i  v  r\  c  r\  1 


Table  2  Calculation  results  of  atmospheric  attenuation 
for  1.315/itn  iodine  laser  (K  “Idem) 


Climate 

tropic 

mid-latitude 

summer 

mid -latitude 
winter  r 

subarctic 

summer 

U.  S.  standary 

text 

otj  km-' 

0.265 

0.187 

0.042 

0.115 

0.075 

0.155 

P,l  km'1 


0.152 

-,-0.123 

0.315 


Climate  conditions  of  the  text  :r**2lt  ,  //“64%. 


extinction  coefficients  are  averaged  up,  the  calculated  j3n=0.137 
agrees  fully  with  the  experimental  value.  However,  the 
experimental  site  was  between  the  urban  and  rural  areas, 
approximately  10km  from  the  city  boundary.  The  aerosol  property 
should  be  between  the  urban  and  the  rural  types .  The  above- 
mentioned  situation  means  that  our  processing  method  is  reliable. 
In  other  words,  the  application  of  atmospheric  transmittance  of 
an  He— Me  laser  to  calculate  the  aerosol  extinction  function  of 
other  laser  wavelength  is  on  a  reliable  basis. 

As  for  the  molecular  absorption  coefficient,  the  average 
value  is  very  close  to  the  calculated  value,  smaller  than  30%. 
There  are  two  possible  reasons;  one  is  the  shorter  transmission 
distance  with  low  instrumental  precision,  thus  causing  certain 
experimental  errors.  From  the  discrepancy  of  data  in  Eq.  (1), 
individual  data  may  have  deviations  of  more  than  30%.  It  is 
possible  to  have  certain  deviations  when  averaging.  Secondly, 
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there  are  some  flaws  in  the  database  of  the  AFGL  spectral  lines. 
Generally  speaking,  the  database  is  reliable  but  there  may  be 
room  for  improvement  due  to  limitations  of  technical  conditions 
a t  that  time  .  In  r 0 c s n t  yesns,  work  in  +■  h -i  o  0 

reported  in  successive  cases.  For  example;  Gao  Junyi  et  al .  [51 

discovered  that  the  ray  intensity  and  1 ine width  given  by  Kilter 
were  on  the  larger  side.  However,  the  results  from  Crossmann  et 
3. 1 .  were  on  the  smaller  side.  If  the  absorption  cross~s8C^1’ An 

v  afi!lecibleP^  *  cmj)  w^ter  vapor  molecules  at  the  1.315- 

micro meter  wavelength  given  by  Bragg  =  ^ ^  1 1  ^  r  1 1  ,, cqa  i -r\ 


ilcu late  the  molecular  absorption  coefficient 


a  +-  r\  V)  1  tr 


0.043km  *  under  the  experimental  conditions.  The  value  is 
reduced  down  to  less  than  half  of  the  AFGL  data.  Further 
discussion  is  necessary  to  determine  which  is  right  and  which  is 
wrong . 

With  this  experiment,  generally  we  can  obtain  the  follcwincr 
results:  in  the  fall  in  the  Hefei  area,  during  the  horizontal 
atmospheric  optical  path  transmission  of  the  iodine  laser,  the 
variation  in  atmospheric  transmittance  is  approximately  between 
70  and  85%.  When  the  ts mp srsturs  is  1  1  ^ -n ^  f v'di  a ^  Tr ^ 

humidity  is  64%,  ths  absorption  coefficient  of  the  atmospheric 
molecules  is  approximately  0.11km”1  for  the  1 . 315— micv*omet£iv' 
iodine  laser.  The  main  factor  affecting  atmospheric 
transmittance  of  the  iodine  laser  is  scatterino  of  the 
atmospheric  aerosol  particles,  and  the  varying  absorption  of 
water  vapor  due  to  variable  water  content  in  the  atmosphere. 


far  end  of  the 
to  measure  only 

the  water  vapor  absorption  coefficient  at  this  wavelength  in  the 
actual  atmosphere . 

To  obtain  more  accurate  transmittance  and  water  vapor 
coefficient  for  an  iodine  laser  in  actual  transmission  in  the 
atmosphere,  we  should  make  experimental  measurements  at  longer 
distances.  In  addition,  we  should  conduct  detailed  studies  in  a 
simulation  laboratory  on  the  relationship  among  temperature, 
pressure,  and  water  vapor  content  at  the  1.315-micrometer 
wavelength  for  the  water  vapor  molecular  coefficient.  Research 
in  this  direction  will  be  pursued  when  conditions  are  suitable. 

The  first  draft  was  received  on  February  26,  1993.  The 
final  revised  draft  was  received  on  Nay  5,  1993. 


Since  the  iodine  1.315-micrometer  laser  is  at  the 


water  vapor  absorption  line,  it  is  very  difficult 
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ABSTRACT  Thermal  blooming  induced  by  aerosol  is  discussed  in  laser  beam  propagation 
though  atmosphere.  We  give  the  effective  absorption  coefficient  ct#//(r).  Because  of  the  limited  con¬ 
ductivity  of  air,  the  increment  of  temperature  of  air  heated  by  aerosol  absorbing  energy  from 
laser  beam  is  imhomogeneous,  which  lead  to  the  dependent  on  time.  The  delay  time  /„ , 

when  the  arf/{i ) arrives  at  the  maximum  value,  is  dependent  on  the  radius  of  aerosol. 

KEY  WORDS  aerosol,  high  power  laser,  thermal  blooming,  effective  absorption  coefficient. 
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In  theoretical  research  on  high-powered  laser  propagation  in 
the  atmosphere,  at  present  research  on  thermal  blooming  is  mainly 
concentrated  on  the  molecular  absorption  effects,  but  neglects 
the  contribution  made  by  aerosol  particles  to  thermal  blooming  in 
the  atmosphere.  If  the  laser  wavelength  at  the  atmospheric 
window  is  properly  selected,  absorption  by  molecules  is  very  low. 
Here  the  effect  due  to  aerosol  particles  should  be  considered. 
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When  there  is  an  absence  of  particle  puncture,  breaking, 
evaporation,  G.  Caledonia  [1]  and  D.  E.  Lencloni  [2]  et  al. 
neglected  the  relationship  of  heat  exchange  between  particles  and 
air.  It  is  assumed  that  all  energy  of  the  aerosol  particles 
absorbed  from  a  laser  beam  is  transferred  to  the  atmosphere  with 
no  time  delay,  the  effective  absorption  coefficient  of  aerosol 
particles  is  the  absorption  coefficient  of  aerosol  particle. 

C.  K.  Chan  [3]  neglected  the  limited  heat  conductivity  of  air;  he 
assumed  that  some  of  the  energy  absorbed  by  aerosol  particles 
from  a  laser  beam  is  used  to  heat  particles,  and  the  remainder  is 
converted  to  heat  and  is  rapidly  diffused  in  the  air.  Thus,  air 
is  uniformly  heated  thereby  inducing  the  effective  absorption 
coefficient  of  aerosol  particles. 

For  the  real  atmosphere,  since  the  coefficient  of  thermal 
conductivity  of  the  atmospheric  molecules  is  smaller,  the  heat 
entering  the  air  is  always  first  stored  in  the  circumference  of 
the  particles,  forming  a  heat  well,  thus  constituting  a  region 
with  a  local  temperature  gradient  prior  to  outward  diffusion. 

Thus,  the  heat  entering  the  air  does  not  immediately  heat  up  all 
the  air.  Research  reported  in  this  paper  aims  at  using  the 
theory  of  the  present-day  continuous-operation  lasers  to  heat 
single  aerosol  particles,  leading  to  thermal  blooming,  and  to  be 
led  into  the  group  of  aerosol  particles  with  consideration  given 
to  heat  exchange  between  particles  and  air,  as  well  as  the  finite 
thermal  conductivity  of  air.  We  study  the  effective  absorption 
coefficient  for  heating  of  group  aerosol  particles  caused  by  a 


27 


q  o  Ci  •y' 


c  a  to 


In  the  second  part  of  the  pacer  ,  the  authors  describe  the 
d  formula  used  by  single  particles  to  heat  the 


surrounding  air.  In  the  third  part,  the  authors  present  the 
models  and  formula  of  the  effective  absorption  coefficient  of  th 
heat  effect  by  group  aerosol 


The  fourth  a  nroconf  c 


the  theoretical  calculation 


II.  Air  Temperature  Distribution  Surrounding  Single  Particles 
after  Heating 

For  convenience  in  this  discussion,  the  authors  neglect 
absorption  by  atmospheric  particles  .  The  temperature  var*'  a  +*  j  r\ m 
distribution  of  air  surrounding  the  spherical  shaped  aerosol 

particles  with  uniform  absorption  can  be  obtained  from  Eg.  (1) 
r  7- r  i 


1  d  T 


=  v:r 


K  0  l 


(1  ) 


Boundary  conditions 


IT 
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cT 


IQ^xa-=  -AnJ  crkj  (  —  +  —  aV  p(cp{  -j—  ) 


T=Tn 


T.tVo  d  TO 


r-  x 


(2) 

(3) 


Initial  condit*’ 


T=  T 


whcin  +■  =  O 


(4) 


Here  k- 

p.c. 


,  K,  ka ,  p2  ,  and  c,  indicate,  respect  ivsly ,  the 


thermal  expansion  rate,  thermal  conductivity,  density,  and 
specific  heat  capacity  of  air.  J  is  Joule's  constant;  and 
the  radius  of  aerosol  particles. 


P  n 


i  to  r\ 


"•  3V>Q  vpcnonfl  TTdl  T7 


the  density  and  specific  heat  capacity  of  the  particles.  I  i« 

o  Q 


the  absorption  coefficient 


the  incident  laser  intensity.  Q;v,c  is 
of  aerosol  particles  to  be  solved  in  the  theory  of  Mie 
scattering.  T  is  temperature  and  Tg  is  the  temperature  before 
the  action  of  air  and  particles.  This  is  under  the  assumption 
that  air  and  particles  are  at  the  same  temperature.  In  Eq.  (2), 
the  authors  have  assumed  that  the  internal  temperature  of  the 
particles  is  uniformly  distributed,  neglecting  the  evaporative 
effect  of  the  particles. 

From  Eqs .  (1)  through  (4),  we  can  obtain  the  approximate 

solution  [3,5]  of  the  aerosol  particles  heated  with  continuous- 
waves  at  a  point  P  (x,  y,  z)  by  the  particles  at  any  points 
' )  in  space . 

2aAT\Aa) 


0(x' ,  y ' , 

T(r,  a  \  r' ,  ()  =  T0  + 

Temperatures  of  particles 


{  1  —  exp  [  — 


t(  a ) 


( 1  -  A-  )]  }exp[  -  /r] du 


(5) 


T  (/)  =  T0+  AT,  (a  ){  1  - exp[  - 


t(  g  ) 


]} 


(6) 


is  the  highest  temperature 


Here  r=|r-r'|.  P=  ,  AT  ^  (a  )  = 

rise  that  can  be  attained  by  aerosol  particles  with  a  as  the 
radius.  ^  is  the  time  constant  of  the  temperature 

increase  of  aerosol  particles.  After  a  single  spherical-shaped 
aerosol  particle  absorbs  laser  beam  energy  and  heats  the  air 
(with  a  as  the  radius),  the  air  temperature  is  given  by  Eq.  (5) 
at  time  t  with  distance  r  from  the  particle  center. 


III.  Effective  Absorption  Coefficient  of  Group  Aerosol  Particle; 
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The  assumption  is  made  that  the  aerosol  particles  a^e 
distributed  with  the  single  dispersion  spectrum;  the  chemical 


r'  ^  w  Ir'  -1- 


es  and 

the  geometric  parameters 

of  all  aerosol 

particles 

s  ame . 

They  are  uniformly  distri 

buted  in  air . 

To 

air  i«. 


are  the  s an 

simplify  the  model,  particle  motion  is  neglected. 

If  the  densi ty  distribution  of  aerosol  particles  in 
uniform,  then  the  probability  of  particle  existence  in  the  volume 
e  1  ement  dx  1  dy  ’  dz  *  near  0  ( x  ’  ,  y  1  ,  z  ’  )  is  d^  1  ^  xr  *  ^  *  fhovofnm 
temperature  rise  of  point  P(x,y,z)  at  a  distance  r  from  0  is 
AT ( r , a ; r 1 , t ) ;  its  probability  if  dx ’ dy ’ dz ’ .  Therefore,  the 
average  temperature  rise  AJ  at  any  instant  can  be  expressed  as 


rci 

L  J 


AT(r,a,t)  =  N\  dx  'dy  ’dz '  2a^TAa) 


(1- 


-exp[ 


r (a)  TT  )]  )  exP[  ~p2]  dp 


(7) 


In  a  unit  volume  near  this  point,  the  heat  required  to 
induce  an  average  temperature  rise  A T  (r,a,t)  of  this  point  at 
time  t  is 

AQ(r,a,t)=JcaptAT(r,a,t )  (8) 

If  we  assume  that  thermal  blooming  of  this  portion  is  uniformly 
absorbed  by  the  atmosphere,  the  effective  absorption  coefficient 


e  *i  rn-i  1 


to  molecular  absorpti on  can  be  used  to  indie  ate  the 


thermal  blooming  inducing  the  average  temperature  rise  of  the 
atmosphere  due  to  the  energy  from  the  laser  beam  abso^b^^  kt?- 


lerosol  particles 
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AQ(0=  nt)*.,fO)dt 


(9) 


Then  the  effective  absorption  coefficient  of  aerosol  particles 
with  r  gspgc  t  to  planar  W3V6S  is 


cap„ 


](t)  dt 


=  f  dx'dy'dz'  2^-  f  -ptepi  -  7 

IU)  J,  Wn  1  T 


-  — (1- A-)]}exp[ 0°) 


IV.  Calculation  Results 

It  is  assumed  that  the  incident  laser  is  planar-wave.  We 
made  computations  on  the  effective  absorption  coefficient  of 
carbon  particles  and  water  droplets  from  a  laser  incidence  to 
single  scattering  spectrum  with  wavelengths  X=2 . Bmicrometers  and 
X=10 . 6micrometers  .  The  particle  concentrations  are  N=10/cm*  and 
particle  radius  a=5micrometers .  The  specific  heat  capacity  (of 
carbon  particles)  cn=0 . 711J/g‘K,  pn=2 . 55g/cm" .  When 

t'  ^ 

X~2  .  Smicrometers  ,  n=l .  8-iO  .  3  .  When  A»=10  .  Bmicrometers  , 
n=2.2-i0.9.  The  specific  heat  capacity  (of  water  droplets), 
c„=4 . 18  J/g’  K ,  pn=l .  Og/cm3 .  When  X=2  .  Smicrometers  ,  n=l .  42-iO  .  013  . 

f  f 

When  A=10 . 6m.icrometers ,  n=l .  13-iO  .  08  .  For  air,  its  heat 
conductivity  ks=2. 54x10  ^J/s.cm.K.  and  its  Joule's  constant 
J=4 . 187 J/cal . 

Figs.  1  and  2  show  the  time  dependent  variation  curves  of 
effective  absorption  coefficient  for  carbon  particles  and  water 
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F.g.l  Effective  absorplion  coefficient  for  Fij.2  Effective  .bsorption  coefficient  .  for 

«‘,b°"  “  3  *"d  ,0-6»"  "■»  droplet  J,!m  „d  ,0  6,; 


coefficient  is  nearly  0  at  0 . 1 microsecond .  With  orolongina  of 

time  ,  the  effective  at>SOrp>^*i/^ri  riAof  f  i  r'l  grit-  c«  I  ^  “|  77-  i  nrivogcoc 

After  a  certain  time  period,  the  coefficient  increases  rapidly, 
until  a  maximum  value  is  reached.  Thereafter,  the  coefficient 
decreases  slowly,  finally  arriving  at  a  steady  state.  Now  this 
is  the  total  absorption  coefficient  of  group  aerosol  particles. 
This  is  a  greater  difference  than  the  results  of  absorption  *Krr 
pure  molecules.  Because  after  the  aerosol  particles  absorb  heat, 
part  of  the  heat  absorbed  will  be  stored  internally  to  form  a 
heat  well,  and  another  part  of  the  heat  heats  the  surrounding 
air.  At  the  beginning,  most  heat  is  stored  i  n  °  ^  <5^  ■hho  ^  o  /o  c  ^  i 
particles;  very  little  of  the  heat  enter0  +-  h  o  0q  ^ -h  f ViQrc 
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is  very  small  variation  of  air 
effective  absorption  coefficient 
the  temperature  rise  of  particle 
the  air  increases.  However,  due 
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to  the  finite  thermal 


conductivity  of  air,  the  heat 
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effect  on  air  variation  increases.  Correspondingly,  the 
effective  absorption  coefficient  begins  slowly  increasing. 
Hereafter,  heat  stored  in  the  inhomogeneous  region  locally 
surrounding  the  particle  also  begins  to  diffuse,  and  the 
effective  absorption  coefficient  increases  rapidly,  to  a  maximum 
value.  Due  to  the  function  of  diffusion  of  stored  heat  and 
thermal  blooming,  and  the  conduction  of  thermal  blooming  from 
particles,  this  maximum  value  of  the  coefficient  is  slightly 
greater  than  the  total  absorption  of  the  group  aerosol  particles 
With  progress  of  thermal  diffusion,  local  inhomogeneity  becomes 
less  and  the  effective  absorption  coefficient  begins  decreasing. 
Finally,  only  thermal  blooming  of  particle  absorption  is 
functioning . 


We  define  that  the  emergence  time  for  the  maximum  value  of 
the  effective  absorption  coefficient  is  the  delay  time  of 
effective  absorption  coefficient  of  aerosol  particles.  To  solve 
for  Eq.  (10) ,  we  can  obtain  the  variation  rate  of  the  effective 
absorption  coefficient  with  time. 


particle  radius  a  /cm 


Fig.  3  Delay  time  tD  versus  particle  radius  a 
( 1  )  water  droplet  ;  (2  )  carbon  dust 


action  of  a  laser  for  carbon  particles  and  water  droplets  with 
different  radii.  From  the  figure  we  can  see  the  following :  with 
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increasing  radius,  the  delay  time  of  maximum  value  of  the 
effective  absorpt  ion  coefficient  increases  Thi0  ic  an co  f 
greater  the  particle  radius,  the  more  the  temperature  rises  to 
reach  eguilibrium.  The  longer  the  characteristic  time,  the  more 
the  heat  stored  surrounding  the  particle,  therefore  the  greater 
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IV.  Conclusions 

In  the  foregoing,  the  authors  studied  the  time-dependent 
variation  process  of  the  effective  absorption  coefficient  for  an 
increase  in  atmospheric  mean  temperature  in  the  process  that 
aerosol  particles  have  the  energy  absorbed  from  laser  beams  to 
transfer  to  air  during  the  interaction  between  aerosol  particles 
and  the  laser  beam.  From  the  calculation  results,  we  can  see  the 
following:  m  the  time  range  below  picoseconds,  the  effective 
absorption  coefficient  of  air  due  to  aerosol  particles  can  be 
neglected.  With  the  time  period  of  picoseconds  to 
submilliseconds,  the  time-dependent  variation  in  the  effective 
absorption  coefficient  of  air  caused  by  aerosol  particles  is  very 
complex.  In  this  time  period,  we  should  adopt  a  method  differing 
from  molecular  absorption  in  processing.  After  submilliseconds, 
the  effective  absorption  coefficient  in  air  (caused  by  aerosol 
particles)  attains  a  steady  state,  not  varying  with  environment. 
This  can  be  applied  directly  so  that  the  overall  absorption 
coefficient  of  aerosol  particles  can  be  used  to  obtain  the 
effective  absorption  coefficient. 
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ANALYTIC  THEORY  OF  SMALL-SCALE  THERMAL 
BLOOMING  INSTABILITIES 

Shi  Jiang jun 

Institute  of  Southwest  Physics, 

P.O.  Box  523-55,  Chengdu  610003 


ABSTRACT  The  linear  analytic  theory  or  small  scale-size  thermal  blooming  instabilities  for  a 
h.gh  energy  laser  propagating  through  a  homogeneous  medium  is  derived  in  paraxial  scalar  wave 
approximation  and  isobaric  supposition.  When  we  perform  Fourier  transforms  in  transversal 
coordinates  and  Laplace  transforms  in  time  and  longitudinal  coordinate,  the  fluctuations  can  be 
obtained  in  analytic  form.  In  the  real  world  the  Fourier  components  of  the  fluctuations  are  written 

w,th  the  propagation  kemal  (or  Green  function)  K„  (z . ,).  after  the  inverting  Laplace  transforms  are 
performed. 

KEY  WORDS  thermal  blooming,  small  scale-size  instability,  propagation  kernel. 


I.  Introduction 

Atmospheric  propagation  of  a  laser  beam,  is  restricted  by 
instabilities  due  to  thermal  blooming.  Small  perturbations  in 
intensity  or  in  phase  of  the  laser  driving  source  will  cause 
modulation  of  high-space  frequency  of  laser  beam  intensity,  thus 
leading  to  variation  in  medium  refractivity  due  to  absorption  and 
other  processes.  These  perturbations  of  the  medium  modulate  the 
laser  beam,  phase,  intensifying  the  rate  of  laser  beam 
perturbation  relative  to  time  and  propagation  distance. 
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Obviously,  this  interaction  between  laser  and  medium  is 
nonlinear .  A  solution  to  this  process  requires  eo nations  from 
fluid  mechanics  on  media  as  x*ell  as  fluctuation  equations  of  the 

Q  jO  O  ~|  ^  Vi  ^  ^  *1  0  ^1  ^  T  ^  V*  O  “J  YY\  T  T  "1  2  f*  ^  ^  -j  ^  ^  ^  q  ^  /^> 

carried  out  by  usinq  data  encodinq  with  a  computer 

However ,  the  field  equations  can  considerably  be  simplified 
m  the  paraxial  scalar  wave  approximation.  In  the  isobaric 
approximation,  the  fluid — mechanical  ecr nations  for  the  medium  are 
very  simple  in  form.  In  this  status,  researchers  can  very  easily 
obtain  the  linear-perturbation  equation  set  relating  to  light- 
intensity  perturbation,  phase  perturbation,  and  medium- 
ref ractivity  perturbation.  In  recent  years  various  *m  "t-  o  o, 
solutions  were  available  for  perturbation  equation  sets,  such  as, 
the  Briqqs 1  analytical  method  of  approachinq  the  solution  [1J 
J.  R .  Morris'  solutions  on  the  standard  Green  function  [2],  and 
T.  J.  Karr's  Laplace  transform  and  inversion  solution  method  [3]. 
Methods  used  in  the  article  are  similar  to  those  used  by  Karr. 
Detailed  derivations  are  presented  in  the  article.  Moreover,  two 
kinds  of  equations  are  given  to  describe  the  propagation  kernel 
K  u(z, t)  for  qro  w  t  h  i  n  p  art  u  r  b  a  1 1  o  n  .  T  h  a  s  a  a  q  u  a^~i/r^T^,c^  •— 1  v*  ci  n  c  0^  ■f 
reference  by  researchers  on  atmospheric  propagation  of  high- 
powarad  lasars 

II.  Theoretical  Fundamentals 

A  laser  field  can  be  described  by  a  complete  sat  of 
Maxwell’s  equations.  The  laser  propaqates  in  the  z -direction ; 

7  Q 


the  propagation  medium  is  nonmagnetic,  with  magnetoconductivity 
equal  to  magnetoconductivity  in  vacuo.  The  electrical  charge  in 
space  and  the  conductive  electric  current  in  the  medium  are  not 
taken  into  consideration.  In  addition,  the  dielectric  constant 
of  the  medium  is  slow  in  time-space  variation.  Therefore,  in  th 
paraxial  scalar  wave  approximation,  the  ligh1"  £  ^ 

6  =  £exp  { i  (u)t-k:) }  (1) 

The  light  field  envelope  E=E (x,y , z , t , )  satisfies  the 
equation 


Vj  E+  2 ik  (  -L  +  )£+  (  J-  -  i  )*’£= 0 

^0 


(2)' 


i., ,,  A,  _uS  the  number  of  laser  waves  and  X  is  the  wavelength  of 
the  laser  waves.  a9  is  the  extinction  coefficient,  which  is  the 
summation  of  the  absorption  coefficient  a,  and  the  scattering 
coefficient  a$;  ag=a3+a£.  Moreover,  a«M0oc/ne;  in  the  equation, 
Mn/  c,  and  n0  are,  respectively,  the  magnetoconductivity, 
electroconductii'i ty ,  and  the  background  refractivity  of  the 
homogeneous  medium,  and  c  is  the  speed  of  light  in  vacuo. 

The  energy  conversion  time  during  which  some  of  the  laser 
beam  energy  is  absorbed  by  the  medium  is  much  smaller  than  the 
propagation  time  of  the  laser  beam;  in  atmospheric  propagation, 
the  energy  conversion  time  is  measured  in  microseconds  and  the 
propagation  time  is  measured  in  seconds.  In  other  words,  the 
* e^ractivity  change  of  laser  energy  due  to  absorption  by  the 
medium  occurs  over  a  long  time  scale.  Therefore,  this  process 
can  be  subject  to  the  isobaric  approximation.  In  the  isobaric 
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medium;  and  D  is  the  mass  diffusion  coefficient  of  the  medium. 
Caused  by  the  unit  laser  intensity,  the  time  variation  rate  f  of 
refractivity  is 

r  =  cba ,  ( 4 ) 


In  the  equation,  c„=  {cnj  cp)  ,/  (ceT) 


DT  =) 


=3  -m  e?  y-\ 

"  ^  v"  ^  ^  Jh' 


here , 


aporoximately  eoual  to  SxlO~^cm"/J  ( 2 n /? P ) ?  is  the  Golstone— Dale 


is  the  isobaric  specific  heat  of  the  medium; 


the  background  temperature  of  the  medium;  and  I  is  the  laser  beai 
intensity 

III.  Linear  Perturbation  Equation 

Assume  In  and  nr  are,  respectively,  the  laser  intensity  and 
refractivity  of  a  medium  in  the  absence  of  perturbation.  Post¬ 
perturbation,  laser  intensity  is  I  and  medium  refractivity  is  n. 
Let  us  also  assume  that  the  pre-perturbation  field  phase  is  $ry 
and  that  it  is  0 1  post-perturbation,  assuming  that  perturbation 
is  very  slight: 


«/„(!  +  /■) 

<pm<p,-<p, 

n  ■"/!,  +  An 


(4.1  ) 

(4.2) 

(4.3) 
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Then  the  light  field  envelope  I  £  I  =  v/  / 


V-\  cz  ur  i  "h  +*  Hi  r\  ^  o 


£=  £0  V  1  +  F  exp  (/ (<p0-<p,)  )*£,cxp  (i>#+  — - /</>,) 


(5) 


Let  us  introduce  the  variation  of  the  relative  ref ractivity 


A  n  1  .  e  , . 

—  ~  or  (  —  -1) 


2  4 

Substitute  Eg*  (3)  into  Eg.  (2) ,  end  substitute  E 
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to  /-t  *1  i  npar  anrwAvi  m  n  h  i  An  \ 

V-*  >>x  X  XiiCui  £-S  J_  X  11.  W  X  f 


dF  1  „j 

(6.1) 

Vi 

(6.2) 

It.  +  v  .  Vu=  -  F+DV*n 

dt  «0 

(7) 

rewritten  as 

^  VJ_  £)  - 0 

(8? 

Eqs .  (7)  and  (8)  are  important  equations  in  the  analytical 

theory  of  thermal  blooming  small-scale  instabilities.  These 
equations  present  the  variation  rule  of  variable-ref ractivity 
propagation  with  propagation  distance  z  and  time  t,  for  a  medium, 
with  the  perturbation  of  laser  intensity. 

To  solve  for  Eqs.  (7)  and  (8),  we  should  present  the  initial 
distribution  fj  (z=0)  of  perturbation  F  (z=Q,t)  for  0,  (z=0,t)  at 
the  boundary  (z=0) . 


IV.  Perturbation  Amount  and  Propagation  Kernel 

Execute  a  Fourier  transform  on  Eqs.  (7)  and  (8) 
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X 


f{x  ,y  )exp(/A  •  r  ^  )dxdyj  (2  r.): 


In  the  equation,  k-ky+kj 
components  of  the  pgv-t n rK: 


,  then  the  equation  of  the  Fouri 
entity  can  be  derived: 


-i  r\  r~\ 


S:F \ 

~jpT  +  a\Ft  =  k-fit 


VFt  —  jQ/r’/i,  +  iv  -  A/j4 


(9) 

(10) 


In  the  equations,  the  transverse— direction  wave  number  k  and 
of  the  perturbation  amount,  and  p‘  are  given  below 

k2  -  k]-r  k\ 

k' 


(11 


2k 

r  =  r  ij  n0 


,j_l  0  0  Q  0  2  p  il  n  C[  i.  FI  l  ^  ^  0  A  p  n  h  "I  /A  TA  3  7A  V"%  'i  T  ^  ^  ^  ^  ^  yi  "!  4"  "1  /A  7A 


f*(/.z  =  0)=Ft(0) 
/•;(/, z=o)=f;(o) 

/h  ( r  ,  t  =  0  )  =  p\  ( r ) 


(12) 


Then ,  carry  out  a  Lad  ace  transform,  on  z~~o  and  t  ~  v  in 
Eqs .  (9)  and  (10),  o  and  v  are,  respectively,  the  conjugate 

variables  of  the  Laplace  transform.  After  manipulation,  the 
express ion  f or mu las  of  the  post-transform  perturbation  amounts 


are  as  follows: 


Vk  = 


~o_  _TK(0)  +  f;(0)) 


*  via :  +  a\)  +  rk:  /i*  /(ffJ  +  a;)  +  'ne 


Fk  = 


k: 


v‘  (<7:+aj)  +  r"A-  v'(ff:ra;)  +  n-: 


rk2 


_  ]  fr-MQ)  +  /V(Q) 


(T-+a; 


(13) 

(14) 
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In  the  equations,  v'  =v -ikv  +  Dkz  •  When  the  convection  velocity  v 
and  diffusion  coefficient  D  are  constant,  the  time  factor 
exp ( -Dk'+ik* v)  is  multiplied  with  the  related  quantities  for 
obtaining  the  inverting  transforms.  We  can  see  that  the  presence 
of  diffusion  has  the  function  of  restraining  perturbation. 

However ,  the  effect  of  the  convection  term  is  to  change  the 
transverse  distribution  of  laser-intensity  perturbation. 

To  carry  out  an  inverting  transform  on  Eqs .  (13)  and  (14), 


1  fH  +“ 


Kl 


introduce  the  following  nronsrfstion  herns! 


K=Kk  (a  ,  v)  = 

L^Lk(a  ,v)s= 

frt  =  Kl(c)m 


Kk,Lk  '  and 


1 


vio'+aD  +  r'k2 

_ <r2+  a\ _ 

vio'  +  aD  +  Tk2 
1 


g,+  aj 


(15) 

(16) 
(17) 


Then  the  equations  for  perturbation  expression  equations  (13)  and 
(14)  can  be  rewritten  as: 


ft r%  ( oFk  (o ) + f/(o  ) ) 


(18) 


Fk  =  k2Ktfil+  ( 1  -  ne^t)  (K”  Fk  (0  )  +  K°k  F'(0 ))  ( 19  ) 

In  the  equations,  K°k'  is  the  Laplace  transform  of  the  . 

Obvi  o  u  s  1  y  / 

~~o _  q  -  (17') 


In  the  above  —  mentioned  equations  r  those  quantities  with  ths 
symbol  ( * )  indicate  a  Laplace  transform  with  respect  to  time. 
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T  h  s  q  vi  ci  n  titles  w  x  t  lx  t  h  £  s  y  jh  x  o  1  *  (  x  5 1  £  5  r  3  t~\  1  3/^0  i-  v  3  n  a  ^  ^  y*  rp 

of  2.  The  quantities  with  superscript  0  indicate  the  initial 
value  of  tins  *  the  independent  quantities  contain inq  0  indicate 
the  boundary  value  of  the  quantities.  For  example,  indicates 
quantity  after  the  z-c  Laplace  transform  of  the  Fourier 
moonents  for  refractivitv  3 x  fv»p  ini  f-i  ai 


+-X 


non 


i  od  , 


Besides  ,  the  quantities  havinq  an  apostrophe  at  the  superscript 


indicate  the  differential  with 


T^gepc/^f  f  n  2 


m x 0  quantities  t*t  -i  t ^ 


indicate  the  differential  of  the  quantity  with  respect  to  t. 
the  definitions  (15)  through  (17)  of  the 


P  -xr  n  c  i  n  a 
^  1 


propaqation  kernel,  we  have  the  foil owinq  relationships: 


cKk  =  K; 

K'^-Vk2  K„  K°' 

-T’k2  KkK\ 

Z4=  -^rj*^(T)jT+i 


(20) 


amonq  other  relationships 


Let  us 


nqrry  All  1 


it  the  invert inq  Laplace  transform  on  z—o  in 
Eq .  (18)  ,  and  the  invert  inq  Laplace  transform  t  —  v  on  Eg  .  (19)  ;  w  e 

can  then  obtain  the  solution  of  perturbations : 

VktZ'O  —  Li'(z<t)n0k(z)—  F  [Kk(z  ,t)Fk(0  ,  t)  +  Kk(z  ,  t)Fk(0  ,  t)]  (21) 

Fk(z  ,  t)  =k2Kt  (z  ,t)n\  ( z)  +  K\'{z)Fk  (0  ,  t)  +  K*k  (z)F'(O  ,  t) 

+  K'(z,t)Fk(0,t)  +  llk(z,t)F'(0,t )  (22) 

W  e  m  u  s  t  n  ote  that  in  the  t  w  o  a  h  o  Tv7  e  m  entioned  equ^t^^n^  n  ti  vu  ^ 
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case  cf  items  containing  two  or  more  pi 


>ducts  ind 


ancn 


p  •h  2 


n  r*  +■  < 


mill  — i  ta  1  "i  V\tt  rr  av  i+'cjvrsc:  PAnfai  ni  nn  f  tja  r>r  ronvci  t\  ^ 

X  XX  K/  J  O  f  Xp/X  -L  WXiUU  VUU  ^U4.iliu^  L  7V  V  X*>  -L  IVtV  i  X  ^  -L  W  X4 

independently  multiplied  by  t,  these  terms  should  be  solved  for 

t*t  ■?  f  Vi  f  V»  q  p  nnTml  ii-f-i  r\  -p  f  In  pn-rcim  prnAv/^i  nrr  +-  r\  f  h  o  i  nncivf  i  nn  T.anl  a^c 

transform.  For  example, 


(r  .  /)F4  (0  ,  t)  -  f Kk  (z  ,  z)Fk  (0  ,  t-x)dt 
Jo 


Kk(z,i)nl(z)  = 


&k(s  *t  )ul  (z-s)ds 


among  others  (not  written  out  here) 

From  the  perturbation  solution  we  know  that  there  are  three 

emir  adc  -f  nr  e  fi  wnl  af -i  nn  -i  p  c?  f  aKi  1  i  f  t;  •  /  1  \  i  nf  dpc-j  f  r^ciV'+-ii-K'V\^>4 — i  r>  ta 

^  w  ui.  x_  w  x_  u  LXiuUx  d  v—  a-  x  x  xix>j  vu  x/  _l.  -i_  -a.  jr  •  \  -x  ;  -i_  a  a  w  x  lx  X  j  ^  i.  vux  vvi 


F,(0,t)  of  the  1 


2  c  or  +-  V' 


qpom-i  eei  An  o  Apr  Ap  (  >-7  =  C\  \  •  /  0  \  nViaop 

CUik>lU  J.  kJ  >J  X  V^JL  U  V  Ux  vC  V  *—*  X/  /  f  \  £~1  /  X 


perturbation  Ft'(0,/)  of  the  laser  transmission  source  (z=0);  and 
(3)  the  initial  perturbation  pk(zr0)  of  refractivity  in  the 
medium.  Only  with  the  existence  of  any  excitation  source  among 
them,  can  the  growth  in  perturbation  (after  interaction  between 
laser  and  medium)  be  given  by  using  the  propagation  kernel 
Kj,  ( z  ,  t )  ,  Lj,(z,t)  and  the  corresponding  post-differential 
quantity.  We  can  see  that  it  is  vital  to  correctly  present  the 
expression  equation  of  the  propagation  kernel  Kj,(z,t). 

V.  Expression  Equation  for  Propagation  Kernel  Kk(z,t).) 

From  the  definition  of  the  propagation  kernel  Kk(o,v)  ,  o-z 
is  conducted  with  the  inverting  Laplace  transform  to  obtain 
K\U.V)  .  Kk(a  ,t)  can  be  derived  with  the  inverting  Laplace 
transform  on  v-t.  There  are  the  following  results: 
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Kt(:  ,v)  -  — —  sin  {aKfi: ) 


(23) 


Kt  (o  ,  t )  =  — r-! — -  exp  ( - ^  - —  ■  ) 

&'  +  a;  cr  +  a\ 


In  the  0QU.3l1.ions 


1  verting  Laplace  t rsns form  of  v~t  is  s^si1" 


ri  /— <  05  ~v~  ~v~  1  czr 


)ut  on  Kk(z,v)  ?  or  the  inverting  Laplace  transform  of  o- 


irried  out  on  P  /  „  ,  ^  +■  h a -n  +■ 

^  ^  Ak  (  O  ,  O  r  - 


Ho  r*  r\y~  v'  o  o  H 


expression  equation  of 


the  propagation  kernel  Kv( z,t)  can  be  obtained.  In  the  inverting 
transform  process,  by  utilizing  the  theorem  of  residues  to  solve 

f  O  JT  t*  1*1  0  ^  0  ^  V*  0  ^  0  Tl  0  f  Vp  11  1  t"  1  —  p  r~\  "1  Q  70  A  *1  TO  t“  1  |  TO  A  4“  "0  AO  TO  C  TiT  O  f“  Vo  4“  Vo  O 


application  of  general  integration 


eg nation  of  the  invertinc 


t  r  p  n  o  f  / 


Ti ,  and  by  utilizing  the  Lagrangian  formula  t 


r\  n  Awn!  q  f  o 
w  s-  wv. 


high-order  differential  of  the  multiplication  product  of  tw< 


functions ,  thus  the  express 


>rously  determined . 


ion  equation  of  K,,  ( z,t)  can 


Beginning  from  Eq .  (23),  first  A',(_-.v)  is  written  in  series 


k -)-£  t-D-'g.--)’-1  _J_  in. 

...  ( 2n+  I  )!  a, 


Then  from  the  above- mentioned  steps  of  the  inverting  transform, 


we  obtain 


M-0=X 


(  — irK.-p-1  j_  y _ «! _ (jkn 

(2n+  1  ) !  at  M%  («-m)!(m!):  a,  1  ; 


nning  with  Eg. 


(24),  Kt(z,t)  can  also  be  obtai 


First,  Kk(a,t)  is  rewritten  in  the  series  form: 


(<f  •  0  =  Z 


( -  i  y(2kr‘att  y  (  1 


n\ 


c2+  a\ 


)'” 


(26) 


There  are  (n+1)  orders  of  pole  points  at  the  virtual  axis  c=±iav 
of  the  a- complex  plane.  Complete  the  inverting  transform  with 
the  above-mentioned  steps,  and  the  result  is 


K  /_  v  (-1  nkrtzy  !  „  (  . 


(26) 


In  the  equation 


(~0"(w  +  m)!  1 

"  m\(n  —  m)\  (2  a**)- 


2  (/)"*"*' 


(26') 


Since  there  are  different  pathways  of  the  inverting 
transform.,  there  are  quite  different  forms  in  K],(z,t)  derived 
from  Eqs .  (25)  and  (26) .  However,  since  the  field  was  begun  fro 

Eq.  (15),  the  results  are  the  same.  If  we  let 

**(*.o-E**w(*.o  ,  we  can  obtain  the  following  from  the  two 

»-0 

equations  mentioned  above: 

K[0)=  —  sin  (atz  ) 

ak 

krmt 

^ri== — —  [a*rcos  (akz) -sin  (o*z)] 

2  (kr't)2  3  3  1 

*J2>  = - - [  -j  sin  (a^z)  -  —  cos  (a*z)  -  (a^sinCa^z)] 


In  the  general  situation,  and  in  the  situation  of  various 
parameters  ( ay ,,  z  ,  t ,  r  *  ,  k )  ,  by  using  the  numerical  dissolving 
method,  we  can  prove  that  the  Ky(z,t)  derived  from  Eq .  (25)  is 

identical  to  the  Kv(z,t)  derived  from  Eq.  (26).  In  Kv(z,t),  the 
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T\  r\  +■  « 


expression  equation  ^ ( 


related  to  tine  is  identical  t  ^ 


the  results  of  the  invertincr  transform  on  ^^e  "*  ^  ^  r^vnr^^rTD  +- -!  nn 

r  1  WJ>  '-'J-X  _i  I_xxv^<^x-  J^/X 

kernel  Eq .  (17)  .  This  explains  that  both  types  of  K|,(z,t) 

expression  equations  are  correct. 
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ANISOPLANATISM  IN  ADAPTIVE  OPTICS: 

EVALUATION  IN  TIME  DOMAIN 

Wang  Yingjian 

Anhui  Institute  of  Optics  and  Fine  Mechanics 
Chinese  Academy  of  Sciences, 

Hefei  230031 

ABSTRACT  In  this  paper,  anisoplanatism  in  adaptive  optics  and  capability  of  phase  derivative 
adaptive  optics  are  analysed  in  time  domain.  The  general  formulas  of  residua!  phase  structure 
function  and  Strehl  ratio  are  obtained,  which  not  only  be  used  to  discuss  the  anisoplanatism.  but 
also  can  be  used  to  analyse  the  other  factors  that  degrade  adaptive  optics  system  performance, 
such  as  the  mean  wind  velocity  and  the  finite  servo  bandwidth  of  adaptive  optics  system.  The  re¬ 
sults  show  that  the  capable  scale  of  phase  derivative  adaptive  optics  are  very  narrow. 

KEY  WORDS  adaptive  optics,  anisoplanatism,  phase  derivative. 


I.  Introduction 

There  has  been  much  discussion  on  anisoplanatism.  in  the 
eddy-current  effect  on  laser  atmospheric  transmission  compensated 
for  with  an  adaptive-optics  system  [1,3].  A  new  method,  the 
phase-gradient  method,  overcoming  anisoplanatism,  was  presented 
in  reference  [3] .  In  the  article,  the  time— domain  analytical 
method  is  applied  to  further  discuss  the  effect  on  compensation 
efficiency  of  an  adapt ive -optics  system  by 

anisoplanatism ,  moan  wind  fisld.  and.  systsiu  bandwidth  Morsovsr 
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T-  In  - 

and 

cAwmcn  c?  a  +-  -? 

■f  V*  AW 

a  unified 

Inn 

^■±x  { 

ent  m e  t  h  o  d  is  (figcjcribed 
domain  ( fre^uency—dom  ai  n  } 
transfer  function  of  the  adaptive-optics  system.  The  article  is 
the  first  to  derive  the  general  equations  of  the  structure 
function  of  the  residual  phase  with  phase  compensation  and 
compensation  under  the  phase  gradient  method.  Next,  the  transfer 
function  of  the  adaptive  optical  system  in  the  n!v 
method  is  derived  in  the  time  domain  of 


‘ace  err  3  r)  -i  on-h 


f  Vic  T7  r 


■  urier  transfc 


+"  ■»  rn  At 


(main .  Finally,  the  results  are  analyzed  and  discussed. 


II.  Strehl  Ratio  and  Phase  Structure  Function 

The  result  of  compensating  for  the  eddy-current  effect  in 
laser  atmospheric  transmission  by  an  adaptive-optics  system 
be  described  with  the  Strehl  ratio  [/!]: 


srp  r*  o  ra 


St~  J dpK(p)exp[  -  j-  D{p)} 


(!  ) 


In  the  ecr nation 


A(p)  = 


f  —  [arccos(  H! -  ) — L£_L  /.  , 

J  *  o,  y-( 

<>  n 


D„ 


■  dpi 


P\^D .) 


(\P\>D0) 


is  the  modulated  transfer  function  of  the  laser  fra 
system ;  Dp  is  the  transmission  aperture;  and  "p 


(2) 


rans mission 
is  the  coordinate 


difference  of  the  deformation  mirror  surf? ( 

OOt  J  C  ciTcf  pm 


50 


D(  p  )  -  Di  ( p  )  +  Dx  ( p  ) 


(3) 


D£  and  are,  respectively,  the  structure  functions  of  the 


fluctuating  phase 


ict  discuss 


nd  oscillation  amplitude.  We  d< 
here  the  effect  on  fluctuating  oscillation  amplitude.  The  phase 
structure  function  is  given  by 

0J(p)=<[A<I>(r,  ,/)-A« >(r2,  t)]}>  (4) 

(r,  r)  is  the  residual  phase  of  the  transmitted  light  beam;  < 

indicates  the  statistical  mean.  In  the  approximate  conditions  o 
geometric  optics,  the  anomaly  phase  generated  by  atmospheric 
transmission  of  a  laser  through  eddy  currents  is  given  by  the 


following  * 


0  — k ^ dzn[  r+  ($z+  v0)  1 1  z  ] 


(5) 


k~2n/X:  the  number  of  laser  waves;  n  is  the  fluctuating 

ref rac tivity ;  h  is  the  transmission  distance;  ^  is  the  angulai 
velocity  of  target  motion  or  the  wind  speed  vector  in  the 
atmosphere;  and  v p;  is  the  mean  wind  velocity  in  the  atmosphere 
As  is  well  known,  there  is  a  time  lag  in  the  adaptive  optical 
system.  Let  us  assume  that  the  time  lag  is  At,  that  the  phase 
relationship  between  the  beacon  light  and  the  main  laser  with 
compensation  of  phase  is  given  by 

Qt(r,t+  At)  =<t»r(r,i) 

®a'®T  indicate,  respectively,  the  phase  of  the  b 
and  the  phase  of  the  main  laser.  Therefore,  during  the  actual 
phase  compensation  there  is  a  residual  phase 


(6) 


o  o  nnn  1  -i  crK  t 

""  a  w 


A<D(r,/+ A/)  =cDr(r,r+ A/)  -4>#(r,  /+  At)  =<t>r(r,  t+  At)  -<Dr(r,r)  (7) 
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Thus  we  know  that  during  the  phase  compensation  th 


o  c  f  y-n  r<  f  i 


function  of  the  residual  ^h 


3  C  C  TO 


T  .  a  +■  7  rr  r 


DJp,  /)  =/d  dz,  dz,  <  [  w  (r, ,  /+  A/,z, )  -n'(r;,  t+At .  z, )  -n  ( r,  ,  t . . 

J  0  ^0 

+  /i(fj ,  f  ,z,)]  x  [n(r, ,  f+ A/  ,Zj)  -  n  (fj ,  /+  A/  ,z}) 

*  with  further  derivation  ,  we  obtain 


(8) 


D,(p,n=k: 


Jo 


|  *[  -  20.(0  ,  z )  +  2Z>.(  |p  | . Z)  +  20.  ( \p  | , 


z) 


-&A\ P  +  P>\'Z)-DA\p-pi\,z)) 


(9) 


In 

the 

eo nation 

pi- 

(0z,+  «•„)  A /  , 

■  o„d 

HIpT+r’)1'1 

*1  o 

the 

st  ructur- 

fun 

c ticn  e 

the 

nor f u vKa 

tion  refract"’  v i t y 

i  s 

the 

s  t  r  n  c  ^ 11  c 

con 

s t an t  of 

the 

nerturba 

t ion  ref rac t iT7i tv 

Since  m  the  general  situation  the  integration  liir.it  of  Az»|p'|,Jr 
can  be  extended  to  (-oc.x)  ,  thus  Eq 


0,(p.  A/ )  =  2.91  fc3 


dzlC;(z])[2\p\il'-\?+pl\,n-\p-pl\>li+2\f,\il']  (10) 


1 


-ft 

V  — 


AC-r)-  [(1  +  x3)l/5+(;cJ)'/,]</;t«2.91.  15 


In  the  above  so nation 
applied . 

If  t h e  p h a s e  g r a d i e n t  c o o e n s a t y  o ^  ^ ^  ^  ^ ^ 


(Dfl(r,f+A()  ~  <£„(')  +  A  b- $,(/)  +  ®Al±*l2_JhiL  A/)  A/  (n) 


Frorr  50°  7  £  )  an/^  /  R  )  t»to 


r<  a  ri  OiKf  a  i  r-> 


#,(**  .  r+ A/ )  *2®r(r,  /)  -<Dr(r,  f-Ar) 


c;o 


(12) 


AO(r,  t+  At )  =k 


dz[  n{r,  t+  At ,  z)  -2n(r,  t , :)  +  n(r,  t  -  A/ ,  r )]  (13) 


By  using  s  s  t  hi  o  cl  siniilsir  to  t h 3. t  rnsnt^  pH  shov 
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a  -n  nKf  a  n 
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MfT)  r-3  Cs-p 

'  i ;  \  ks  ^  i  i 


"i  B  n  3  S  0 


OAp.At)  f 


nr  the  rSSiduHl  pH a eg 


D,  (p',At)=  2.91k 


cJ  |*rfr,C;( 
Jq 


)[6i^is/j-4ip+/Tli'j-4i/r-^; 


1 5/ 3 


+  I  ip  +  2p~t  |5/3+  I  "p  —  2p^  |5M+  8 1  pi  I5 1  -  2|  2/T,  Is  j] 


(14) 


In  t li 0  foil owincj f  wo  3  p  p>  1  y  this  ti  —  b— cio'osin  Pouir i on  - 1 r s n. s f or^ 
method  (as  described  in  reference  [2])  for  processing,  that  is, 


4»r(r,/)  =  J*rf/<t>r(r,Oexp(  —  271 rift] 


(15) 


f  is  the  conjugate  variable  of  the  Fourier  transform  for  time  t. 
<j )s(r,/)  can  be  used  to  describe  the  following  relationship: 

< t>e(r,f)=H(f)<t>T(rJ )  (16) 

H ( f )  is  the  Fourier  transform  of  the  time-domain  transfer 
function  in  the  adaptive-optics  system.  With  processing  similar 
to  that  mentioned  above,  the  general  equation  of  the  residual- 

phase  structure  function  [2]  can  be  obtained: 

... 


D,(p  )  =  2.91  k- 


dfdtdz,  |  1  -//(/)  F  exp  [-2nifi] 


•  c;U)[I|p+p,r+ 


(17) 


In  the  phase  compensation  situation,  from  Eq.  (6)  we  know 

tf(/)«cxp[-2*»yAf]  08) 

Let  us  substitute  Eq .  (18)  into  Eq.  (17)  and  we  can  obtain  Eq . 
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(10) .  With  respect  to  compensation  under  the  phase  gradient 
method,  from  Eq.  (12)  we  can  obtain 

<DS  (r,/)=2tf(/  )<t>T(r,f)-If(f  )Or(r.  /)  =  //,  (r,/) 
that  is,  during  compensation  under  the  phase  gradient,  method,  the 
Fourier  transform  of  the  time- domain  transfer  function  in  the 
adaptive-optics  system  is 

H,(J)  =  H(f  ){2-H(f))  (19) 

Thus,  let  us  substitute  Eqs .  (18)  and  (19)  into  Eq .  (17),  and  we 

can  very  easily  obtain  results  similar  to  that  in  Eq .  (14) . 

III.  Discussion 

From  the  foregoing  analysis  we  know  that  Eq .  (17)  is  the 

general  equation  describing  the  structure  function  of  the 
residual  phase  for  the  time— domain  error  in  the  adaptive— optics 
system.  More  detailed  discussion  [2]  was  conducted  on  the  effect 
of  anisoplanatism  in  the  situation  of  response  bandwidth  and 
phase  compensation  in  the  adaptive— optics  system.  Here  we 
discuss  mainly  the  difference  of  the  effect  by  the  atmospheric 
mean  wind  field  and  anisoplanatism.  in  phase  compensation  and 
compensation  under  the  phase  gradient  method.  Furthermore,  we 
can  know  the  adaptation  range  in  the  phase  gradient  method. 


3.1.  Effect  of  Mean  Wind  Field 


When 
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reduce  the  Strehl  ratio,  we  first  analyze  two  extreme  cases. 

When  lpj/r0|«Z)0/r0  from  Eqs .  (1),  (10),  and  (14),  we  obtain 

the  Strehl  ratios  St  and  St,  of  laser  atmospheric  transmission  of 
phase  compensation  and  compensation  under  the  phase  gradient 
method,  respectively,  as: 


St  =  exp[  -  6.881  —  |5/ ^ 
ro 

5/,  =  exp[  -  5.676  I  |5/5  ] 


(20) 


In  the  equations,  r e  is  the  transverse-direction  coherent  length. 
From  Eq.  (20),  then  we  can  see  that  compensation  under  the  phase 
gradient  method  is  better  than  phase  compensation.  We  define  it 
as  follows:  when  the  Strehl  ratio  is  reduced  to  1/e, 

IPil  is  the  isoplanatic  length.  Then  in  the  situation  of 
phase  compensation  and  compensation  under  the  phase  gradient 
method,  the  isoplanatic  lengths  pn  and  pn1  are,  respectively: 

p0  =  0.3144ro 

p0l*=  0.3528  r#  (21  } 

We  know  that  the  isoplanatic  length  with  compensation  under  the 
phase  gradient  method  is  1.122  times  the  isoplanatic  length 
during  phase  compensation. 

When  I p,/ r„|» D0/ ra  ,  that  is  the  case  of  entirely  different 

optical  paths  for  the  beacon  light  and  the  main  laser,  then  we 

have  , 

D.(p)a 2x6.88|-£-  |5'j 

-  •  r*  ■/■*':  ■ 

A  (p) -6x6.88|~  pj 
rB 


However,  the  phase  structure  function  without  phase  compensation 
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is  6 . 8 6  Ip?/-,!’'1  ,  then  we  can  see  that  the  phase  compensation  is 
not  worse  than  in  the  absence  of  compensation,  and  compensation 
under  the  phase  gradient  method  is  even  worse  than  phase 
compensation.  In  other  words,  applying  an  unrelated  phase  for 
compensation  in  the  present  transmission  situation  is  even  worse 
than  the  absence  of  compensation. 

For  the  general  situation,  from.  Eqs .  (1),  (10),  and  (14)  we 

can  obtain 


nD\ 


Si,= 


nDl  j 


</p7:(p)cxp{  -  3 .44  [  2 1  —  |,M  - 1  —  }  - 1  — — |5,J+  2|  — 15;,1 } 

r«  ro  ro  ro 

* 

(J'pk('p)ex p{  -3.44 [6|—  |>/J-4|  — +  [5/J  —  4 1  —  Z£i  p  > 


(23) 


+  1 1± ~£j  |5/i+  |  £zlfh  |»'J+  1>65|  £l  |»/J]  } 
ro  ro  ‘ 


Fig.  1  shows  the  comparison  of  the  Strehl  ratio  for 
atmospheric  transmission  of  a  laser  in  the  presence  of  eddy 
current  during  phase  compensation  and  with  compensation  under  the 
phase  gradient  method  in  different  situations  of  jp^/rol  ‘ 
a,  b,  c  are  the  results  of  phase  compensation  ;  d,  e,  f  are  the 
results  of  compensation  under  the  phase  gradient  method. 
a‘  c!  .  I  pj  r0 1  0.1,  b,  e  ,  I  pj  ra  |  -  0.3 1  ;  c,  f :  I  pj rQ  |  *  1,0.  .  From  Fig  .  1  ,  we  can 

see  that  the  results  of  I  pj  r0  |«  DJ  r0  and  |  ~pjr0 1»  A)/ro  are  consistent 
with  the  results  in  the  above-mentioned  analyses.  For  different 
\pjr9\  ,  there  is  a  critical  D^/r^.  When  D^/r^  is  lower  than  the 
critical  value,  in  the  case  of  compensation  under  the  phase 
gradient  method  the  decrease  in  the  Strehl  ratio  with  increase  in 
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Fig.  1  Strehl  ratio  v*  DJ  r,  at  different  I  p’j  rf\  Fig.  2  Strehl  ratio  vr  DJ  r0  at  different  0/  0„ 

with  phase  and  phase  derivative  compensation  with  phase  and  phase  derivative  compensation 


Dg/rg  is  more  rapid  than  that  under  phase  compensation,  and  it  is 
smaller  than  the  Strehl  ratio  in  phase  compensation.  When  Dg/rn 
is  greater  than  the  critical  value,  in  compensation  under  the 
phase  gradient  method  the  Strehl  ratio  more  rapidly  approaches 
the  limiting  value;  this  is  the  result  shown  in  Eq.  (20).  Then 
compensation  under  the  phase  gradient  is  better  than  phase 
compensation.  Moreover,  with  increase  in  lp,/r0|  ,  the  critical 

value  of  Dg/rg  is  higher.  The  physical  reason  is  as  follows:  on 
the  one  hand,  when  !/T,/r0|  is  a  constant,  with  increase  in 
Dg/rg,  the  residual  phase  is  greater  in  phase  compensation  and  in 
compensation  under  the  phase  gradient  method;  thus,  the  Strehl 
ratio  decreases  with  increase  in  Dr,/rn .  On  the  other  hand,  when 
Dg/rg  is  smaller,  with  rg  held  constant,  | "pj r# I  is  greater,  and 

lp,/Z)#|  is  also  greater.  The  phase  relationship  between  the 
beacon  light  and  the  main  laser  is  more  deteriorated,  thus  not 
only  is  the  residual  phase  difference  greater  under  phase 
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P  i*l  3,  S  S  3*.  Cl  S  n  t  TH  0  od  3,  S  1.  •  1  2  2  to-  TP  0£  +■  Vi  T  o  at>  1  a  p  a  t-  n  r  a  t~i  ri  "1  o  nn/^  cv 

phase  compensation.  However,  when  0 » Dj  h  >  there  is  the  similar 
result  as  in  Eq.  (20  )  .  In  o  t  bi  s  u  words  ,  th.s  phsss  co^pGnsst ion  is 
worss  then  in  tins  s b s s n c s  of  cornpons s t ion  snd  compensst ion  Linds r 
the  phass  gradisnt  msfcbiod  is  svsn  worss  than  ohass  con^psnsation 
In  the  aeneral  situation  f mm  T?n-  / 1  ^  m  r»  ^  =>  r->  ri  m/i\  tjq 

calculate  the  Strehl  ratio  in  phass  compensation  and  in 
compensation  under  the  phase  gradisnt  method  for  different  O/©^ 
values.  The  results  are  shown  in  Fig.  2a,  b,  c,  as  the  results 


in  phase  compensation;  d,  e,  and  f  are  the  results  in 
compensation  under  the  phase  gradient  method.  a,  d  :0/0o  =  0.5; 

b.  e  :  0/0„  =  l.o  ;  c,  f:0/0o  =  5.O  •  The  calculation  condition  is  that 
there  is  the  horizontal  homogeneous  atmosphere.  With  respect  to 
an  inhomogeneous  atmosphere,  the  results  do  not  intrinsically 
differ.  From  Fig.  2,  the  fundamental  result  is  similar  to  the 
effect  of  the  main  wind  field. 

From  the  foregoing  analysis,  only  in  the  case  when  there  are 
smaller  \  p\l  rq\  and  ©/© g,  and  D^/rj  is  greater,  the  compensation 
under  the  phase  gradient  method  then  be  better  than  phase 
compensation.  From  Figs.  1  and  2,  we  see  that  when 

•P'/'J  =  0.1  and  0.31,  as  the  phase  gradient  method  is  better 

than  the  phase  method  that  the  critical  values  of  Dj/r^  are, 
respectively,  20  and  55.  When  ©/©g=0 . 5  and  1.0,  the  critical 
values  Djj/rg  are,  respectively,  30  and  70.  In  the  case  of  even 
smaller  ©/©g  and  I  <pj  rB  |  ,  there  is  very  slight  differences  in 

the  effect  on  compensation  under  the  phase  gradient  method  and  on 
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in  the 


phase  compensation.  With  the  greater  |p,/r0|  and  ©/©,., 
effective  range  the  phase  gradient  method  does  not  have 
superiority.  Even  when  D^/r-  is  much  greater  than  |pj/r0|  ,  or 

Dr,/r-  is  much  greater  than  ©/©,-,  the  isoplanatic  length  and  the 
isoplanatic  angle  in  compensation  under  the  phase  gradient  method 
are  only  1.122  times  the  isoplanatic  length  and  the  isoplanatic 
angle  in  phase  compensation. 

4.  Conclusions 

By  using  the  time  domain  analytic  method,  the  author  and  his 
colleagues  obtained  the  general  equations,  Eqs .  (1)  and  (17),  of 

the  residual  phase  structure  function  and  the  Strehl  ratio  in  the 
effect  of  an  adaptive-optics  system  during  anisoplanatism ,  mean 
wind  field,  and  the  system,  finite  bandwidth.  The  results  are 
consistent  with  those  in  reference  [2] .  The  effect  on  mean  wind 
field  and  anisoplanatism  in  phase  compensation  and  in 
compensation  under  the  phase  gradient  method  are  emphasized  in 
the  analysis.  As  revealed  in  the  results,  when  \Jf,+  v  |A t«Dli 
the  isoplanatic  region  compensated  under  the  phase  gradient 
method  is  only  1.122  times  the  isoplanatic  region  compensation  in 
phase  compensation.  However,  when  I  Oh +  v J A/  ,  the  phase 

compensation  is  worse  than  the  absence  of  compensation,  and  the 
compensation  under  the  phase  gradient  method  is  even  worse  than 
phase  compensation.  For  general  cases,  even  as  the  Strehl  ratio 
in  compensation  under  the  phase  gradient  method  is  several  times 
greater  than  that  of  phase  compensation,  but  then  the  Strehl 
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ratio  is  considerably  lower  than  the  actual  requirement.  We  can 
see  that  the  compensation  of  the  eddy-current  effect  in  laser 
atmospheric  transmission  is  not  very  effective  under  the  phase 
gradient  method.  Even  to  a  certain  extent,  compensation  under 
the  phase  gradient  method  is  better  than  phase  compensation,  but 
the  effective  range  is  very  narrow,  and  the  extent  of  improvement 
is  also  very  limited.  Moreover,  there  will  be  greater  complexity 
in  the  adaptive-optics  system,  and  higher  computation  volume  in 
the  wavefront  processing. 

The  author  expresses  his  gratitude  to  researcher  Song 
Zhengfang  for  valuable  discussions  during  the  writing  of  this 
paper . 

The  first  draft  was  received  on  February  15,  1993;  the  final 
revised  draft  was  received  for  publication  on  June  5,  1993. 
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ABSTRACT  The  verification  of  testing  power  of  ground  -  based  high  power  laser  is  an  impor¬ 
tant  PhysiCAl  problem  in  arms  control.  In  this  paper  we  propose 'a  scheme  to  measure  the  scat¬ 
tering  light  of  laser  beam  in  atmosphere  by  multi -detector.  The  optimal  detective  parameters  are 
calculated  and  several  algorisms  to  evaluate  the  laser' power  are  studied. 

KEY  WORDS  verification  of  laser  power,  arms  control,  scattering  by  aerosol. 


I.  Introduction 

The  firing  level  of  high-powered  laser  weapons  is  directly 
related  to  the  power  of  laser  devices.  To  prevent  further 
development  of  high-powered  laser  weaponry  from  becoming  a  threat 
to  space  targets  and  to  inhibit  the  arms  race  in  space, 
scientists  proposed  the  concept  of  limiting  the  development  of 
laser  weapons  by  limiting  the  power  at  which  lasers  can  be  tested 
in  the  atmosphere  [1] .  This  requires  an  appropriate  method  of 
evaluating  laser  power  in  testing.  Although  it  is  convenient  and 
effective  to  directly  measure  laser  power  in  laboratories,  such 
evaluation  schemes  may  be  difficult  to  accept  by  those  under 
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inspection.  Therefore,  a  method  should  he  found  to  estimate 
laser  power  from  a  certain  distance.  T.  H.  Braid  [2,3]  proposed 
a  scheme  to  estimate  laser  power  by  measuring  the  aerosol 
scattering  light  in  the  atmosphere.  Under  the  assumption  that 
scattering  is  isotropic,  Braid  et  al .  also  proved  that  the 
scattering  light  for  a  laser  beam  capable  of  attacking  targets  in 
outer  space  is  very  intense.  There  is  still  a  very  high  signal- 
to-noise  ratio  detectable  beyond  distances  of  1km.  Therefore, 
this  evaluation  scheme  is  feasible. 

Actually,  the  aerosol  scattering  of  laser  beams  in  the 
infrared  band  is  anisotropic;  the  difference  of  light  intensity 
at  different  scattering  angles  can  be  assigned  in  three  orders  of 
magnitude.  Since  the  direction  of  the  test  laser  beam  can  be 
randomly  varied  within  a  certain  range,  the  scattering  angle  also 
changes  from  laser  beam  to  sensor.  Therefore,  the  power  of 
scattered  light  received  by  the  sensor  is  also  variable  so  there 
will  be  great  error  in  estimating  the  laser  power  based  on  this 
approach.  To  reduce  the  estimation  error  engendered  by  uncertain 
laser  beam  direction  to  meet  the  requirements  of  evaluation,  the 
authors  proposed  a  scheme  for  detecting  scattered  light  with 
multiple  sensors. 

During  target  practice  tests  with  lasers,  the  target  can 
also  be  used  to  evaluate  laser  power  from  the  diffused  reflected 
laser  light.  The  evaluation  method  discussed  in  the  paper  does 
not  include  this  case. 
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II.  Fundamental  Principle 

Usually,  the  atmospheric  transmission  and  the  comprehensive 
demonstration  tests  of  laser  weapons  are  selected  at  sites  with 
good  meteorological  conditions,  conducted  on  sunny  days  with 
higher  visibility.  Here,  laser  beam  transmission  over  thousands 
of  meters  will  not  suffer  from,  obvious  attenuation.  Moreover, 
the  function  of  multiple  scattering  can  also  be  neglected.  The 
sensors  used  to  monitor  laser  scattered  light  is  in  the  order  of 
magnitude  of  approximately,  in  distance  from,  the  laser  beam,  in 
the  range  of  kilometers,  much  greater  than  the  laser  beam,  range. 
Under  these  conditions  the  power  of  laser  scattered  light 
received  at  the  sensors  is: 

Po=P 

In  the  equation,  P  is  the  laser  beam  power.  3(r,0)  is  scattering 
coefficient  of  the  volumetric  angle  at  point  r  with  scattering 
angle  0a.  a  is  the  incident  angle  of  the  scattered  light  at  the 
sensor.  S  is  the  sensor  area.  1  is  the  distance  along  the  light 
beam  from  any  point  r  to  the  light  source.  D  is  the  distance 
from  point  r  to  the  sensor. 

Within  the  range  of  thousands  of  meters  near  the  ground,  the 
aerosol  particle  size  distribution  is  almost  constant,  with 
concentration  gradually  decreasing  with  height  h.  Therefore,  the 
scattering  coefficient  of  the  volumetric  angle  can  be  expressed 
as  : 

P(r,e)=p(9)cx  p(--£~)  (2) 


J 


P(r  ,6)cosaSdl/  D2 


(1  ) 
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In  the  equation  3 (©)  is  the  scattering  coefficient  of  the 
volumetric  angle  for  the  aerosol  on  ground  level;  and  H„  is  the 
datum  of  the  aerosol  [4,5] 

//,« 1.4km  (3) 

Hence, 

'Wsj#<f’)exp(- -L  )^-dl  (4) 

If  the  aerosol  particle  size  distribution  in  the  laser  beam 
direction  as  well  as  the  sensor  position  and  direction  are  not 
changed,  then  power  Pj  of  the  scattered  light  detected  is 
proportional  to  laser  power  P.  By  using  a  laser  with  known  power 
as  a  scale,  laser  power  P  can  be  estimated  based  on  scattered 
light.  The  uncertainty  of  laser  beam  direction  will  cause  errors 
under  this  estimation  method.  The  estimated  related  error  is 
correlated  with  the  shape  of  the  scattering  coefficient  3 (©theta) 
of  the  volumetric  angle.  Fig.  1  shows  three  typical  scattering- 
angle  curves  on  sunny  days  [6] .  From  Fig.  1,  we  can  see  that 
this  scattering  is  anisotropic. 

III.  Evaluation  Scheme 

Along  the  circular  circumference  at  a  distance  Xr  from  the 
light  source  to  be  evaluated,  several  sensors  are  placed  at  equal 
spacing  along  the  circumference.  And  each  sensor  is  facing  the 
direction  of  the  light  source.  The  elevation  angle  of  the  sensor 
is  y.  Its  longitudinal-direction  opening  angle  is  <o;  its 
transverse-direction  opening  angle  is  1805.  So,  the  angle  of 
f ield-of-view  is  u  x  1803. 
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Fig.  1  Angular  scattering  functions  of  water 
haze  with  Young  particle  -  size  distribution 
function  (p^^O.OSpm  .  pmt *  5.00pm  . 

power  exponent  Pt  ) 

’  Curves  I  :  3  ,  x-  10.6pm 

2  ‘  P,m  3  ,  2  •  1 . 1  5pm 

3  :  P  =  4,  i»  10.6pm 


Let  us  assume  that  the  position  of  the  laser  light  source  is 
at  the  origin  of  coordinates .  The  connected  line  between  the 
light  source  and  any  sensor  i  is  the  x-axis.  Then  the 
coordinates  of  sensor  i  are  (xgf0,0) .  Assume  that  the  pole 
angles  of  the  laser  beam  are  (0ff,  #t)  r  then  the  coordinates  of 
any  point  on  the  laser  beam  are 

(/sinfl.cosg,,  Is\n0msin<p„ ,  IcosOJ  .  The  height  of  point  r  is 

h  =  /  costfm 

After  calculations,  the  distance  from  point  r  to  the  sensor 

i  is 

D=  (x20+  I2  -2.x0lsinOmcos<pm)'n-  (6) 

The  scattered  angle  from  point  r  scattering  to  sensor  i  is 

COS0=  *0  Sing,  sing.-/ 

(.x„+/:-2.r0sing,coscp,  )''2 
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The  angle  of  incidence  of  the  scattered  light  onto  the 


sensor  i  is 


COSO!  = 


cosy(-t(,-/sin0.lcos</>.,)  +  /  cosO^siny 
(xl+  i: -  2x0lsinOmcos(pmy  - 


(8) 


Since  there  is  a  certain  limitation  at  the  longitudinal- 
direction  for  the  f ield-of -view  angle  for  the  sensor,  points  of 
too  high  and  too  low  on  the  laser  beam  are  situated  outside  of 
the  f ield-of-view  angle,  therefore  light  passing  through  these 
points  cannot  reach  the  sensor.  By  calculation,  we  know  that  the 
lowest  point  capable  of  reaching  the  sensor  is 

/  =  — _ £n _ _ 

mm  - -  /  q  k 

c°s0„+  sind^cos^ctg  (y-  ) 

the  highest  point  capable  of  reaching  the  sensor  is 


oc  ,  ^cos0„,+  sin0mcos<pmctg(y  +  —  )^0 


cos0„+sin0wcos<p„,ctg(y+  ~  )  and  others 


(10) 


Let  us  substitute  Eqs .  (3),  (5),  (6),  (7),  (8),  (9),  and 

(10)  in  Eq .  (4)  and  we  can  calculate  that  the  scattered  light 

power,  which  can  be  received  by  the  sensor  i  is 

r  ^*** 

P(0)e\p( --JT-  )  ^r-  dl  (II) 

J I  np 


A  calculation  was  made  by  the  authors  by  using  an  example  of 
a  1km  distance  between  the  sensor  and  the  laser  source ,  taking 

(12) 


.t0=  I  km 
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V.  Power  Estimation  Method 

Based  on  Eq.  (11) ,  we  can  determine  the  scattered  light 
power  measured  for  each  sensor.  When  the  laser  beam  direction  ^  «=■ 
upward ,  the  scattered  light  cower  measured  bv  each  sen°or  ^  c 
same.  When  the  laser  Voa  arn  cl-1  rcmf  i  nn  cViancTqc  rfhnDvol  1  17  cr\G3il*-i  mrr 

there  are  different  scattered  light  power  values  received  at  each 
sensor.  The  authors  tried  three  algorithms  to  average  the 
scattered  light  power  received  by  each  sensor.  Then  the  average 

?d  to  estimate  the  laser  beam 


scattered  light  power 
(a)  Ari thmeti 


n  c  1 1  o  o  r 


TN/OT.TCi'r 


mean 


PD=  -  (IPo.,  ) 


(13) 


In  the  equation,  n  is  the  number  of  sensors.  is  derived  from 

Eq .  (11). 

{ b )  P o s t - s e 1 e c t i on  average 

Remove  the  maximum  and  minimum  power  of  the  i  -i  rr>r(- 

ust  measured  ,  the  arithmetic  cgsc  is  f  o 1 1  n ^  -f >" -m  -g c 
(n-2)  data: 

vr  1 


ncs  rcmai  m  nA 


po=  -—2  f  Ipo./-max(?0i()  -min(P0  ,)]  ,  3 


(14) 


(c)  Weighted  average 

A  certain  order  of  power  of  PD  <  is  used  as  the  weighting 


factor  theref c 
average : 


TO  f~*  O 


n  also  be  called  the  oower — b  = 


(15) 


In  the  equation  w  is  a  parameter  that  can  be  used  to  calculate  an 


optimal  value. 

After  not  a  long  period  (less  than  lh)  before  and  after  the 
laser  tests  in  the  atmo  sphere  for  the  evaluation  a  laser  beam 
with  known  power  P*  is  vertically  emitted  f ro^  the  laser  light 
source  in  order  to  scale  the  system.  Now,  the  power  ^  .* 
at  each  sensor  is  the  same.  The  following  eguati 


£  C  6 1 V  6  ( 

on  can  be  used 


to  e s t r m ate  1 ^s^r 


.  "  j 

n»'  n  *  r 


Pmp0- 


0.0 


Tn  -t-V-io  Q/-Tii  a  f  "i  ca  ri 


,  n  i  c  rfi  non  H  t  r  T?rrcr 

P  D  -*  ^  J. 


(13) .  (14 


(16) 

iie;i 


The  estimated  1  as 
actual  lase 


p  V'  TA  rA  T.7  O  V  D  * 


nr  nr\ 


t  be  definitely  egual  to  the 


The  es t i™at°d  v-  o  ~i  a  n 


TTjzs  q  v  r  n  r  n  e? 


B- 


p'-p 


(17) 


(  Q  )  avo  nnf  i-  "h  o  a-mc 


The  estimated  errors  from,  different  laser  beam  directions 

T3  A7  nh  armi  nrr  Vi  o  i  va  r<  +•  i  nn  v  n  n  ol 

*-*  JT  '»"■**  V-'*-  HXJ.  UiVU*  Y  i  i  ^ 

igle  of  the  sensor,  the  error  v a 1 

-ie  a  ■F-iirA/^-h-i/'A-n  /-a  -F 


n  a  c  t, 


rji  1  1 


a  I  o 


change .  In  other 


words  € 


o. « v. . y 


.7/  -r^^v  ,*f. 


(18) 


We  assume  that  the  included  angle  between  the  laser  beam  and 
the  vertical  direction  does  not  exceed  60°  because  the  distance 
of  a  laser  beam  passing  through  the  atmosphere  is  0F  exceeds  60*. 
This  makes  it  disadvantageous  for  a  laser  beam  to  pass  through 
the  atmosphere .  The  authors  reguest  that  the  estimated  errors  € 
of  all  possible  laser  beam  directions  and  powers  should  be  as 

c :  q 


Fig.  2  Error  varies  as  the  angle  of  laser  beam  <pm 

ssnsors  3.ncl  *y — £  3 .  3 v  Vis  c sn  S00  that  the  in  0  3.  r.  p  o  is  si  a 

rn  inimum  wiisn.  ibis  field-~of — v  i  e  w  sriy  1  s  Ci)  aporoaches  0.  Ac  tusll  v  , 
the  f ield-of -view  angle  cannot  be  equal  to  0.  Only  with  a 
smaller  f ield— of —view  angle,  to  be  given,  such  as  (0=0. 001"  or 


0=1 


wg d n  err or  still  a^^roach  the  ^in*^  ir*um 


T7  2  I  11CC 


r\  v'  f  Vi  c  r«  h  anrrej 


Iff  av  -f-  ^  c<  -1  n  rr  i-  Vi  d  Vi  r-i  nr  Ci  r\  -F  "1  p  c*  Ci  v*  t.T2  T7C  1  onaf  l-i 

li.4.  W  A.  V*  V*  V  V I-  c_  J.  J.  vj  UiJ  ■ 

xn  the  aerosol  particle  size  distribution  with  climate,  the  shape 

nrvp  will  r> 


of  the  scattering^ angle 
direction  and 

arlTran  no  c?  r\  f  K  occ  r*  a 


"hanged.  However,  the 


anrtl  cs  a1 

uuy  c 


T71  Ci  T.T  r\ 


•F  +-  K 


Ci  O  Ci  T>  C?  V"  3  V“  O  OOt-  T 


Ci  anfhrxV'c  o 


HIT 


ar\nl  n  o  c? 

^  ^  _L.  J-  >_  V^. 


f  the  field  of 
nnot  be  changed.  Th€ 

■hjxB  angle  scattering  curve  X  to  attain  three  sets  of  detection 
parameters  in  order  to  calculate  the  mean  errors  of  the 

onaffer-i  nrt  —  an/tl  c  nn  vttoc  TT  an^  TTT  t  r  p  nl  n\  r>  +-  -}  n  nr  +-  In  v*  ci  ci  nnucr 

*—>  ■■ —  VC  V  v  w  -k_  J.  ii  V-J  ci  ii  x  v  v  cli.  v  v  v  ci  i  x  cx  j.  x.  x-  •  i_i  jr  d  cX  v  ^  v_Li.iv)  tux  c.  u  ^  tv  c _ _ 

estimation,  methods  and  three  scattering-angle  curves,  different 
mean  errors  and  maximum  errors  are  obtained  for  different 
sensors.  These  values  are  listed  in  Table  1. 


VI.  Conclusions 

Prom  Table  1,  we  can  see  that  the  mean  error  can  be  s— aller 
than  20%  and  the  maximum  error  can  be  smaller  than  50%  due  to 
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